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1.0 INTRODUCTION 

The condensation of vapor on the inside surface of a tube is an important 
heat transfer mechanism in many thermal devices. An example of such a device 
is the heat pipe. This is a liquid-vapor two-phase system, often in the form 
of a closed tube, which transfers heat by evaporation and condensation at 
opposing ends of the tube. The heat pipe has been used to advantage in many 
aerospace and satellite systems, and is now being recognized as a potentially 
attractive heat transfer device for terrestrial applications. The intentional 
or accidental introduction of noncondensable gases into heat pipes signifi- 
cantly changes the condensation heat transfer characteristics of the heat pipe 
condenser. This unique behavior has been widely used for temperature control 
in gas-loaded heat pipes. This report describes experiments and analyses done 
to further establish the characteristics of these gas zones in the limit as 
the gas becomes increasingly soluble in the working fluid. 

Early analyses of the gas-controlled heat pipe assumed a very simple model 
for position and effect of noncondensable gas within the condenser; that is, 
the noncondensable gas was assumed to be in the form of a plug with a sharply- 
defined vapor-gas interface. In time it was found that the flat-front model 
did not accurately represent many systems. It was observed that axial heat 
transfer via the condenser wall caused the sharp front to become more diffuse 
and that the boundary was also elongated due to vapor diffusion into the stag- 
nant gas zone. The culmination of one-dimensional models including wall con- 
duction was the numerical program GASPIPE^^^ which allowed the user flexibility 
in modeling condenser heat transfer effects and yielded better agreement be- 
tween design calculations and experimental observation. 
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In this study a more exact one-dimensional condensation heat transfer 
model for insoluble gases has been developed and compared with experimental 
data. Modifications to this model to accommodate soluble gas behavior have 
also been accomplished, and the effects on gas front behavior demonstrated. 

In Section 2, Analysis, analytical models for condensation heat transfer 
are documented. In Section 3, Experimental Apparatus and Technique, a novel 
optical method used for measuring gas concentration profiles is outlined. 
Section 3 also presents experimental data. Section 4 is an interpretation of 
these data and Section 5 is the summary and conclusion. 
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2.0 ANALYSIS 

In this section models are developed for condensation heat transfer and 
dissolved gas transport in a gas-loaded heat pipe. The influence on condensa- 
tion of an insoluble noncondensable gas is developed in Section 2.1, while 
various effects of a soluble gas are explored in Section 2.2. The intention 
of these sections is to construct an analytical framework that can be used for 
comparison with experimental heat transfer data presented in Section 3.8, 
Experimental Results. 

2. 1 One-dimensional Model for Condensation in the Presence of an Insoluble 

Noncondensable Gas' 

To determine the effect of gas solubility on condensation heat transfer 
in a heat pipe, it is worthwhile to develop an analytical model for condensa- 
tion in the presence of an insoluble gas so that a bench mark or reference 
model is available for comparison with data. An additional incentive for 
developing such a model evolved during execution of this study--the soluble 
gas models, in addition to being complex, were numerically unstable. These 
problems are discussed in more detail in Section 2.2, Soluble Gas Modeling. 

The coupled heat transfer and mass transfer processes in the condenser 
have been described by two one-dimensional conservation equations accounting 
for heat and mass transfer in the vapor phase and heat transfer along and 
through the heat pipe wall. Both axial and radial pressure gradients are 
neglected, and composition and temperature are assumed unifonn across the 
entire vapor cross section. The boundary conditions used are similar to those 
of Marcus, but certain simplifying assumptions have not been made. Specif- 
ically, the first and second spatial derivatives for vapor and wall tempera- 
ture are not assumed equal, allowing accurate representation of high thermal 
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resistance systems such as the glass heat pipe used in experimental phases of 
this study. The modeling described in this section is similar to modeling by 
Saaski,^^^ but various improvements have been incorporated which Increase 
solution accuracy and reduce computational time. 


2.1.1 Vapor Phase Mass Transfer 

The intent of this modeling is to define heat transfer in the stagnant 
gas zone shown in Figure 2.1 and existing between the axial positions 
z<zsz . A cylindrical heat pipe cross section is assumed, but the 
equations to be developed are not geometry-dependent. However, uniform wall 
and vapor cross-sectional areas are necessary. An immiscible gas/fluid combi- 
nation is also assumed. If axial and radial pressure gradients and thermal 
conduction are negligible in the vapor phase, and the heat pipe wall is uni- 
form in temperature along any circumferential line, then vapor migration into 
the gas zone can be modeled as diffusive flow into a stagnant gas layer. Con- 
servation of working fluid vapor (species a) at any cross section results in 
the coupled heat and mass transfer equation 


■ar 




v«fg 


(2-1) 


where denotes molar flux and T and T represent the vapor and wall tempera- 

V W 

tures respectively. is the molar heat of vaporization and the vapor cross- 
sectional area is A^. This expression is based on the vapor phase being sepa- 
rated from the wall by a liquid film of uniform conductance per unit length, 

G^. This constant thickness film is attributable to a uniform thickness capil- 
lary wick upon the wall. In the case of reflux operation without a wick, the 
film thickness is a slowly varying function of axial position, but can often 
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FIGURE 2.1. REPRESENTATIVE THERMAL PROFILES FOR A GAS-LOADED HEAT PIPE 

be approximated by a constant value. The surface of this film which is in 
contact with the vapor phase is at the vapor phase temperature; any interfacial 
resistances are assumed small. For an insoluble gas, there is no net gas move- 
ment and (|)^ = 0 (b denotes the noncondensable gas). The molar flux of species 
a is given as 

( 2 - 2 ) 


A.. H,. 


Tnrx^ 


so that 


V ■ ' - X, } 


(2-3) 
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The ultimate goal is to obtain a differential equation which has as 
the primary dependent variable. With this objective in mind, the mole frac- 
tion X. in a constant pressure system is expressible as 

d 


(2-4) 


W V 


(2-5) 


where is a function of the vapor temperature T^, and is the total heat 
pipe internal pressure. The concentration-diffusivity product, can be 

expressed as a function of T^ through the Chapman-Enskog formula for dilute 
gases 


CDab ' 


(2-6) 


where B is a very slowly varying function of temperature, and T^ is expressed 
in absolute temperature. For practical purposes, B can be evaluated at a mean 
temperature and assumed to be a constant. 

When the expressions above are substituted into equation (2-3), a second- 
order nonlinear differential equation in T^ results 
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(2-7) 


The numerical technique used to solve this equation is discussed in Section 
2.1.3. 
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2.1.2 Axial Wall Conduction 

The heat pipe wall is assumed uniform in cross section and to have a 
temperature dependence on the z-coordinate only. The wall is coupled to the 
vapor phase through a film or layer with per-unit-length thermal conductance 
and to the heat sink through a similar conductance, G^. In Figure 2.2, 
this series of thermal conductances is shown as both a nodal network and as 
related to the condenser-heat sink assembly. A steady-state heat balance on 
an infinitesimal wall slice leads to the differential equation used in the 
experimental phase of this study: 



where T is the heat sink temperature and A and 
s w 

sectional area and thermal conductivity. 


are the wall cross 


2.1.3 Numerical Solution Method - One-dimensional Problem 

In the last two sections, second-order differential equations for T^ and 
T have been derived. Under steady-state operating conditions, both equations 

W 

must be satisfied simultaneously at all axial positions. 

By formation of a nodal network as shown in Figure 2.2, it has been pos- 
sible to solve the coupled equations (2-7) and (2-8) by converting the equa- 
tions to finite difference form and using successive over-relaxation (SOR) to 
iteratively approach consistent profiles for T and T . However, the approach 

V W 

has required long computational times and alternate methods for solving the 
equations were investigated. All failed but one. Most notable among the 
failures was the so-called "shooting technique" in which initial guesses for 
the temperatures and their derivatives at z * z^ are adjusted until the 


4 


f 



SYSTEM COMPONENT 


VAPOR TEMPERATURE 
FLUID FILM CONDUCTANCE 
WALL TEMPERATURE 

INTER FACIAL CONDUCTANCE 
HEAT SINK 


DESIGNATION 


EQUIVALENT FINITE DIFFERENCE NODAL NETWORK 


T„(VAPOR) 


FLUID 

FILM 

(Of) 


GAS GAP 



T (CONDENSER 
WALL) 


T, (HEAT SINK) 


CROSS SECTION OF 
HEAT REJECTION ASSEMBLY 


FIGURE 2 . 2 . FINITE DIFFERENCE MODEL FOR GAS-GAP CONDENSER/HEAT SINK ASSEMBLY 



r 

I 




I 


9 


solutions at z = obtained by integrating the differential equations across 
the zone, satisfy a zero derivative boundary condition at z^^j. However, the 
solution for would invariably either rise or fall with great speed after a 
number of apparently stable increments in z. The nonlinear equation (2-7) has 
many of the mannerisms of a "stiff" equation. 

Successive over-relaxation is comnonly employed as a point-iterative 
technique; that is, one point at a time is adjusted in value. A block relaxa- 
tion technique was defined which was stable and improved convergence rates. 
Taken individually, the finite difference forms of (2-7) and (2-8) are tri- 
diagonal matrices for T^^ and T^^. , if temperatures on the other nodal line are 
fixed. That is, in the nodal network of Figure 2.2, if the wall temperatures 
are held fixed at one iterative step, then the system of linear equations de- 
fining the axial vapor temperature profile is tridiagonal. A tridiagonal 

(31 

matrix can be solved exactly using the Thomas algorithm.' ' This exact set of 
vapor temporal ijres is then held constant, and a new line of wall temperatures 
is calculated using the Thomas algorithm. The new wall temperatures are held 
constant and a vapor profile is calculated again as in the first step, complet- 
ing the iterative cycle. All theoretical profiles were generated using this 
block iterative technique. 

Prior to generating difference equations, the temperatures and axial co- 
ordinate system were nondimensionalized as follows: 


(2-9) 

where T^^^ is the adiabatic vapor temperature. The vapor pressure f(T^) was 
found by fitting experimental vapor pressure data to the following form: 
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• Unf = (2-10) 

. The vapor temperature is in absolute degrees. The finite difference form 

of equation (2-7) can be configured into a tridiagonal matrix with each row 
given by 


For intermediate points 
by the following set of 


^i \,i * ^i \,i+l 


excluding the boundaries, 
equations: 


d. (2-11) 

the constants are defined 


“1 






-(2 + F ,) 



(2-12) 

(2-13) 

(2-14) 

(2-15) 

(2-16) 

(2-17) 

(2-18) 


(2-19) 
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The finite difference equations for the wall were similarly derived and 
will not be presented. Boundary conditions used in modeling were 


T* = 0.995 

T* = 1/(1 + G,/G^) 




+ + 

" < ^0 

+ + 

z = 



0.75 


+ 8.5 

0 


( 2 - 20 ) 

( 2 - 21 ) 

( 2 - 22 ) 

(2-23) 


Standard methods were used to incorporate these boundary conditions into 
the finite difference equation sets; for example, the derivative boundary 
conditions were expressed as central divided differences. 

The dimensionless vapor temperature was calculated over the range 
zt < z^ < zt + 8.5. At a distance of 8.5 dimensionless units from the gas 
front, the gas concentration is essentially constant and calculations over a 
longer zone are not necessary and reduce solution accuracy in the area around 
Zq for a fixed number of increments. The axial gas zone was divided into 114 
increments. Wall temperatures were calculated over a wider range extending 
into the nonblocked condenser zone. This is done because the heat pipe wall 
in the gas zone acts as an extended fin, draining heat away from the warmer 
areas and depressing wall temperature beyond the actual gas front. 

A listing of the computer program used for modeling purposes is given in 
Appendix C. The code is in BASIC language and was programmed on a PDP-11 
minicomputer. Representative profiles are tabulated in Appendix B, and a 
comparison with experimental data is given in Figure 3.10. Iteration of 
temperature profiles was terminated when the largest chc.ige in wall 


12 


temperature per iteration was less than 0.001 in dimensionless units. This 
typically required 40 iterations--at the end of these iterations, the vapor 
temperature was often very slightly less than the wall temperature for 
» z^. With a larger number of iterations this physically unrealistic 
situation is corrected, but the difference in profile shapes is quite negli- 
gible, and was not considered worth the additional computer run-time. Forty 
iterations required about 5 to 6 minutes of CPU time on the PDP-11. 

2.2 Soluble Gas Modeling 

An initial attempt to model condensation in a tube containing soluble 

gas was done using a von Karman-Pohl hausen integral technique. By assuming 

characteristic equations for the concentration and velocity profiles in both 

vapor and liquid phases, it was possible to generate coupled differential 

equations describing the two-dimensional concentration and flow fields for 

incompressible flow. The approach used was similar to that of Somogyi and 
(41 

Yen.' ' The analysis developed will not be presented here in detail because 
the authors were unsuccessful in obtaining stable numerical solutions. 

A shooting technique was used from the interface at z = z^. A stable 
solution set was obtained until the differential integration approached the 
fully stagnant gas plug zone. Under these conditions, the fluid film on the 
heat pipe wall is saturated with a large amount of noncondensable gas, and 
the gas-saturated liquid moves toward the evaporator into regions of 
lower vapor-phase gas concentration. This supersaturation creates a large 
gas flux out of the fluid surface and, in fact, the integral techniques indi- 
cated total vapor phase fluxes that were directed radially inward. When this 
occurred, the solution became unstable and was not correctable with smaller 
steps in z. Because of these large radially inward gas fluxes, it was not 
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possible to obtain a zero derivative for concentration and zero value for the 
axial and radial velocity components at z = z^. The integral method was 
abandoned and the hybrid 10/2D method described in the following section was 
developed. 

2.2.1 Hybrid Soluble Gas Model 

For tubes on the order of 1 cm diameter and vapor pressures of 1 

atmosphere (typical heat pipe conditions), vapor phase diffusivity can be 

much higher than liquid phase diffusivity. Vapor phase diffusivities range 

2 

typically from 0.1 to 1.0 cm /sec, while liquid phase diffusivities of 1 to 
-5 2 

10 X 10 cm /sec are typical for gases in common working fluids. 

Therefore, if radial concentration gradients are neglected in the vapor 
phase but are accounted for in the liquid phase, it is possible to extend the 
analysis of Section 2.1 to include soluble gas effects for a wide range of 
physically relevant heat pipe designs. 

If interfacial resistance is neglected at the vapor/liquid interface 
as in Section 2.1, then the concentration of gas at the liquid phase inter- 
face is in the equilibrium with the gas concentration in the vapor phase, and 

^gt “ ® ^gv liquid/vapor interface) (2-24) 

where C , is the liquid phase gas concentration, is the vapor phase gas 
gx gv 

concentration, and a is the Ostwald coefficient. A mass balance on the 
working fluid vapor yields the same expression as equation (2-1). However, 
the noncondensable gas is no longer stationary, and the equation (2-2) for 
4), must be modified to 
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'as 


*a*b - '^ab^a/^ 

T-*a ■ 


(2-25) 


and the conservation equation for coupled heat transfer and vapor transfer Is 






(2-26) 


To obtain the flux (j)^^ in (2-25), transport of noncondensable gas within 
the surface liquid layer must be detailed. In the experiments described in 
Section 3, the heat pipe was wickless and vertically oriented with the 
evaporator at the bottom. The surface layer therefore is in th-^ classic 
falling film category as investigated by Nusselt and others. If it is 
assumed that the dissolved gas has no effect on physical properties of the 
liquid and that diffusion of dissolved gas in the vertical direction is negli- 
gible compared to convective transport, then conservation of gas in the film 
is governed by 


V(y,z) 




(2-27) 


where V(y,z) is the film velocity and is the diffusion coefficient of the 
gas in the liquid. The coordinate y is measured radially inward with respect 
to the heat pipe inner wall surface. 

Under conditions of laminar flow, negligible shear at the liquid-vapor 
interface, and negligible vapor density compared to liquid density, the 
velocity profile is 



(2-28) 



The fluid density and viscosity are and respectively. The fluid 
film thickness is 6. At the wall, V(0) is zero, while at the vapor/liquid 
interface, V(l) is a maximum. 

Boundary conditions on are 





(2-30) 



At 2 = the dissolved gas concentration has reached a maximum value in 

equilibrium with the fully stagnant gas concentration, 

The liquid film thickness 6 changes with axial position because of 

condensation. By using the velocity distribution (2-28) and assuming a 

linear temperature profile across the film, film growth is obtained by 

integrating over the heat transferred from z = Zj^ to z = z^, where z^ z 

z . The film thickness 6 is given by 
m 



(2-33) 


where t. is the heat pipe wall thickness, 
w 

For coupling to the finite difference equations it was necessary to 
include the wall thermal resistance between the fluid film and the node at 
the center of the wall. For metal-walled heat pipes this radial wall 
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temperature drop can be neglected, but It is of first-order importance in a 
glass heat pipe because of a much lower thermal conductivity. 

For numerical solution of the conservation equation (2-27), a dimension- 
less gas concentration and radial measure were defined. 



The gradually-changing film width 6 must be accounted for in the con- 
servation equation 


V(y*,z) 



d« 

dz 







(2-35) 


The conservation equation (2-35) with boundary conditions 2-30 through 
2-32 was solved by implicit finite difference methods similar to those 
described in Section 2.1. The gas phase concentrations used in boundary 
conditions 2-30 were obtained from a finite difference solution of the vapor 
phase equation (2-26). This solution was also obtained in a manner entirely 
analogous to the solution method in Section 2.1 except for inclusion of the 
new term in ()>., i .e. , 

d 

X 

new term = (1)^^ (2-36) 


The axially-varying flux was obtained from the solution to the liquid- 
phase diffusion problem. At each axial position, a mass balance including 
both the vapor and liquid phases shows that the flux in the vapor phase 
must be balanced by a gas flux of opposite sign within the liquid phase. 
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The liquid-phase gas flux is obtained by numerically integrating the 

concentration profile, 


I = aTrD^6 


z-z 


1 




( 2 - 37 ) 


where is the heat pipe inner diameter. 

The coupled problem was iteratively solved by alternately holding the 
solution fixed in one phase and solving the finite difference equation in the 
other phase. At the completion of each cycle, would be updated. When 

the concentration profiles in each phase remained constant with iteration, 
the coupled problem was considered solved. 
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3.0 EXPERIMENTAL APPARATUS AND TECHNIQUE 

3. 1 Alternate Methods for Measuring Vapor Phase Gas Concentration 

To compare experimental data with analytical models for condensation in 
the presence of a soluble gas, an approach is desired which indicates 
accurately both vapor phase composition and spatial variations in composition. 
Many previous investigators have inferred internal vapor phase dynamics from 
variations in heat pipe wall temperature, or from probe temperatures in the 
vapor phase. The first method has the obvious deficiency of not being a 
direct measurement and, in fact, it is easy to show that measured axial 
temperature variations are primarily related to axial wall conduction for 
most metal-walled heat pipes. It is very difficult to infer from external 
measurements an accurate picture of vapor phase mass transfer in the transi- 
tion zone between the freely condensing and fully stagnant zones. 

It is also difficult to accurately measure the condition of the vapor 
phase directly. Any probe inserted into the vapor generally has such a large 
thermal mass that profiles can be severely distorted. Thermocouple techniques 
or thermistors are generally difficult to integrate into an operating heat 
pipe, not only affecting the accuracy of temperature measurement through their 
own mass, but also fundamentally changing the flow distribution. 

If a small diameter thermocouple well is placed along the axis of a 
cylindrical heat pipe, boundary conditions for the Navier-Stokes equations 
are modified as: 

BEFORE PROBE: at r = 0, U * U^, 

AFTER PROBE: at r = 0, U = 0, 


3r 


= 0 


1^-0 

3r 
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where U represents the axially-directed velocity component. This is a very 
fundamental change in velocity distribution and, therefore, the measurement 
technique has significantly influenced the system being measured. 

3.2 Vapor Composition by Optical Spectroscopy 

As discussed, spatial variations in noncondensable gas concentration are 
difficult to measure within an operating heat pipe without obtaining a dis- 
torted picture of vapor composition. This difficulty has been overcome to a 
great extent by using optical spectroscopy. 

The basic concept is shown in Figure 3.1. A working fluid/gas combina- 
tion is selected in which the noncondensable gas is either strongly absorbent 
or optically transparent at a wavelength where the working fluid vapor is of 
opposite characteristic. Using a monochromatic light source, :uch as a laser 
or spectrophotometer, the outgoing beam from the source is shaped by a narrow 
slit perpendicular to the heat pipe axis. The heat pipe envelope is glass, 
and light transmitted through the heat pipe impinges upon a sensitive photo- 
cell coupled to a pi coammeter. The heat pipe/condenser assembly is mounted 
on a vertical drive so that the heat pipe can be scanned by the beam. 

Absorption of light by a gas is given by the familiar Lambert-Beer Law 

I = I„ EXP(-BC(,z^j,) (3-1) 

where is the incident light intensity, p is an absorption coefficient, and 
Zab the distance of light travel through the absorbing medium. is the 
effective, geometrically-weighted average concentration along the path at any 
given axial position. At constant wavelength a photo-diode has a linear 
response to light intensity, so that I and can also be regarded as photo- 
detector current. Therefore, if p and z - are known, concentration of the 
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FIGURE 3.1. BASIC TECHNIQUE USED FOR MEASURING CONCENTRATION PROFILES 
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absorbing species, C^, can be found. 

^b “ 6^ (t) 

The constant a now is a cell constant and assumes uniform geometry along the 
heat pipe length. 

This technique allows very accurate measurement of gas concentration 
versus axial position with minimum disturbance to the heat pipe vapor. 

Spatial variation along the heat pipe axis can be defined to any degree of 
accuracy by decreasing slit width, although practical limits arise in obtain- 
ing adequate light intensity to drive the photodetector. 

The technique also has some limitations. It is obvious that the 
technique is principally useful in defining the average concentration as a 
function of axial position; the beam in effect integrates any radial varia- 
tions. It is also not advisable to use any wick upon the heat pipe inner 
diameter, as it will optically scatter light. Therefore, this technique is 
most useful for studying condensation in the presence of noncondensable gases 
upon smooth vertical surfaces. Such surfaces are of considerable importance 
in two-phase heat transfer equipment; for example, many terrestrial appli- 
cations use heat pipes in vertical reflux, in exact agreement with the method 
optimum for testing. In addition, both vapor phase and condensate phase 
conditions are very similar to conditions in zero-g, where a wall wick and 
surface tension pumping replace the falling condensate film defined by viscous 
shear and gravity, 

3.3 Apparatus Description 

Table 3.1 summarizes the experimental system. 
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TABLE 3.1. EXPERIMENTAL TEST VEHICLE 
CHARACTERISTICS 


Heat Sink 
Type 

Length 

Diameter 

Optical slot width 
Optical slot height 
Heat pipe centering method 

Conductance of centering 
screws 

Heat Pipe 

Envelope type 
Wall thickness 
Nominal length 
Evaporator length 
Adiabatic length 
Working fluid 
Gas additive 

Optical Assembly 

Operating wavelength 

Effective beam size ac 
heat pipe 

Nominal indicator current 
with no absorbing gas 

Background with light 
source shielded 

Resolution 


Gas-gap chiller block with optical 
slot; Aluminum 

15.24 cm 

1.4618 ±0.0013 cm (I.D.) 

0.318 cm 
13.20 cm 

Four 6-32 nylon screws at each end 
of sink 

0.012 W/K for each set of four 
screws 


1.395 cm O.D. Pyrex glass tube 
0.12 cm 

50.0 cm 

16.0 cm 
18.8 cm 
CCI4 

Chlorine 


3550 A 

0. 10 X 0.318 cm 

4.03 (10"^^) ampere 

^0.02 (10"^^) ampere 

±0.025 (10'^^) ampere on 10"^® scale 


i. 
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Figures 3.2a and 3.2b present perspective views of the glass heat pipe 
and aluminum heat sink. The heat sink, which was in the form of a split-shell, 
had an inner diameter of 1.416*0.0013 cm, and was 15.24 cm long. Four cool- 
ing channels of 0.63 cm diameter maintained the condenser isothermal. The 
optical transmission slits were 0.318 cm wide and 13.2 cm long. After the 
heat sink was machined, it was given a dark brown chromate conversion coating 
having a high absorptivity in the blue portion of the spectrum. The outside 
was then painted flat black. In Figure 3.2b, the light beam is visible in 
the optical slit. 

Before operation, the glass heat pipe of 1.345 cm O.D. was inserted into 
the inner hole, and centrally located using a feeler gauge and 8 radially 
mounted nylon 6-32 screws; four screws are at the top of the condenser and 
4 at the bottom, spaced at 90-degree intervals. These screws served to hold 
the heat pipe in the condenser with a minimum conductive loss, and establish 
a uniform gas-gap circumferentially. 

The heat pipe was 50 cm in overall length. The evaporator was 16 cm 
long. To maintain the evaporator and adiabatic sections as thermally isolated 
as possible, an aluminum walled Dewar, continuously pumped with a nitrogen 
trap and fore pump, surrounded these sections, and was sealed to the heat pipe 
at the condenser end with a Buna-N 0-ring. Electrical feed-throughs supplied 
current to the evaporator resistance heater through the opposite end of the 
Dewar. On the external surface of the Dewar, a uniformly-space resistance 
heater maintained the exterior of the Dewar within 1°C of the heat pipe adia- 
batic temperature. This was necessary to reduce thermal losses via the 0-ring 
contact at the condenser end, and via a mounting plate at the opposite end 
which centrally located the heat pipe in the Dewar. Tiipse elaborate pre- 
cautions were necessary since the heat pipe was operated as much as 75°C above 
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FIGURE 3.2a. PERSPECTIVE VIEW OF TEST SYSTEM USED IN MEASUREMENTS OF AXIAL 
NONCONDENSABLE GAS CONCENTRATION DISTRIBUTIONS. BY USING AN 
OPTICALLY ABSORBING NONCONDENSABLE GAS. A PRECISE BEAM OF 
MONOCHROMATIC LIGHT CAN DEFINE CONCENTRATION PROFILES. 
VARIATIONS IN TRANSMITTED LIGHT INTENSITY ARE DIRECTLY RE- 
LATED TO CHANGES IN GAS CONCENTRATION. 


FIGURE 3.2b. 


FRONTAL VIEW OF CONDENSER AND HEAT SINK SECTIONS. 
TRANSMITTED LIGHT BEAM IS VISIBLE NEAR THE BOTTOM 
OF THE HEAT SINK GROOVE. 
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laboratory ambient. The condensjgr was capable of moving vertically via a 
friction drive mechanism so that any portion of the transmission slit could 
be intercepted by the detector beam. 

3.4 Gas and Fluid Selection 

Based upon a high absorption in the blue spectrum, chlorine was selected 
as the noncondensable gas. A typical absorption spectrum for CI 2 gas at 
1 atmosphere and 20°C is shown in Figure 3.3. An operating wavelength of 

O 

3550 A was selected for all optical tests. This wavelength is very close to 
the absorption peak of Figure 3.3. As a working fluid, carbon tetrachloride, 
CCl^, was chosen on the basis of compatibility with chlorine, no significant 
blue spectrum absorption for thin fluid films, and a high affinity for taking 
chlorine into liquid solution. The latter characteristic is essential to 
define the effect of solubility on gas zone behavior. The solubility of Cl 2 
in CCl^ and R-11 was experimentally measured using an optical absorption 
technique. Results of these experiments are shown in Figure 3.4. 

3.5 Calibration; Optical System 

It was experimentally established that the logarithm of optical absorp- 
tion was not linearly dependent upon gas concentration. Experimental 
absorption measurements with known concentrations of chlorine in a 1.4 cm 
diameter cell established the following relationship between gas content in 
standard cm^ per cc and tn (Iq/I) at A = 3550 A and 20.8°C. 

Cg = 0.39815 in (I^/I) + 0.066075 tn(^)<1.35 (3-3) 




/ 


A 








FIGURE 34. THE SOLUBILITY OF CHLORINE IN Ca^ AND R-11 (CFCI 3 ) 

C = 0.29267 - 0.30541 4n (l^/l) * 0.49546 (in 

/'o\ 

I 1.35 

The experimental data basis for these expressions is shown in Figure 3.5. 

At a higher temperature T, the concentration is 

Cg(T) * 1 - [o. 0005441 «.n(I^/I) (T-20.8)] •Cg(20.8°C) (3-5) 

The light source employed at all times was a 500-watt slide projector 
operated from a constant-voltage transformer. The wavelength selected for 

o 

operation, 3550 A, was produced with an EG&G monochromator. Model 585-22. 
Beam thickness at the heat pipe was chosen at 0.1 cm as a compromise between 
spatial resolution and photodetector resolution. The radiometer detector 
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£n (l^/l) 


FIGURE 3.5. EXPERIMENTAL CORRELATION OF GAS CONCENTRATION WITH THE 

LOGARITHM OF THE OPTICAL ABSORPTION RATIO FOR A CYLINDRICAL 
CELL WITH IDENTICAL PHYSICAL DIMENSIONS TO THE TEST HEAT PIPE. 



30 


head and indicator units were also manufactured by EG&G, and were Models 
580-22A and 580-llA, respectively. 

3.6 Thermal Measurements 

Heat sink and heat pipe temperatures were measured using Type K thermo- 
couples. Adiabatic temperature was measured with a thermocouple centered in 
the adiabatic section on the heat pipe exterior wall. All temperatures were 
measured with 0.1°C resolution using a Doric DS-350 digital thermocouple 
indicator. 

3.7 Data Reduction 

At completion of a test series, data included heat pipe vapor core 
temperatures and heat sink temperatures, heat input data, and photometer 
current readings as a function of axial position in the heat pipe. 

Heat losses or gains through the Dewar surrounding the evaporator and 
adiabatic sections were negligible. However, the nylon locator screws at 
the condenser entrance form a direct heat leak to the heat sink, and this 
heat loss was subtracted from the gross heat input to obtain net heat de- 
livered to the condenser. Each set of nylon screws had a combined conduc- 
tance of 0.012 W/°K. 

Gas concentrations were calculated using Equations (3-3) through (3-5), 
with the sink temperature used as the characteristic temperature in equation 
(3-5). The reference photocurrent was based on the peak current immedi- 
ately in front of the gas front. For most data runs, this current was about 
4.03 X 10"^® ampere. Total gas content within the gas zone was calculated 
by numerically integrating the concentration profile with Simpson's Rule. 

Gas concentration beyond the farthest measureable condenser position (about 

13.7 cm) was assumed equal to the concentration at 13.7 cm. 
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3.8 Experimental Results 

Concentration profiles from experimental sequences using Cl 2 /CCl^ at 
sink temperatures of 14°C and 40°C are shown In Figures 3.6 and 3.7, 
respectively. The axial position scale Is with reference to the condenser 
end nearest the evaporator; e.g., the blind end of the condenser Is dis- 
placed 16 cm axially from the "front end" of the condenser. The continuous 
curves, which represent theoretical predictions, are discussed In Section 4. 
The raw and converted data are both given In complete tabular form In 
Appendix A. 

Vapor temperature response to heat Input Is displayed In Figure 3.8 
for both sink temperatures. Figure 3.9 presents the gas zone molar charge 
based on Integration of the experimental concentration profiles In Figures 
3.6 and 3.7. The molar gas load Is shown versus stagnant zone gas concen- 
tration to determine If gas Is going Into solution In the stagnant zone. 
Figure 3.10 presents two models for predicting noncondensable concen- 
tration distributions compared with data for the 2.2 W high temperature 
CCl^ run. 

Analysis of the data presented here Is discussed In the following 
section. Data Interpretation. 
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FIGURE 3.6. A COMPARISON OF EXPERIMENTAL CONCENTRATION PROFILES WITH THE THEORETICAL ONE- 
DIMENSIONAL MODEL ASSUMING NO SOLUBILITY EFFECTS. GAS CONCENTRATION PROFILES 
WERE CALCULATED FROM OPTICAL ABSORPTION DATA. SINK TEMPERATURE - lA.^C 


' 'O' 



NONCONOENSAM.E CONCENTRATION (t«d. 




VAPOR-SINK TEMPERATURE DIFFERENCE 


70 


60 


50 


40 


30 


20 



0 2 4 6 8 10 12 14 16 18 20 

HEAT TRANSFER RATE (WATTS) 


FIGURE 3.8. A COMPARISON OF CALCULATED TEMPERATURE HEAT TRANSFER CHARACTERISTICS WITH EXPERIMENTAL 
DATA. CALCULATED VALUES ARE BASED ON THE KNOWN VAPOR-SINK TEMPERATURE DIFFERENCE AND 
THE NUMERICALLY INTEGRATED TOTAL GAS CONTENT. THIS INFORMATION IS INPUT TO THE ONE- 
DIMENSIONAL CONDENSATION HEAT TRANSFER MODEL DESCRIBED IN ANALYSIS SECTIONS. AND A GAS 
PROFILE AND HEAT TRANSFER RATE ARE THEN CALCULATED. 
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FIGURE 3.9. TOTAL AMOUNT OF NONCONDENSABLE GAS CONTAINED WITHIN STAGNANT GAS 
ZONE AS A FUNCTION OF STAGNANT ZONE GAS CONCENTRATION. GAS CONTENT 
OBTAINED BY NUMERICALLY INTEGRATING EXPERIMENTALLY-DETERMINED 
CONCENTRATION PROFILES. NOTE THAT THE VAPOR PHASE APPEARS TO LOSE 
CHLORINE GAS AS THE GAS CONCENTRATION INCREASES. 
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4.0 DATA INTERPRETATION 

4.1 Comparison of Experimental Results With Theory 

The full complement of data profiles was simulated with the one- 
dimensional insoluble gas model and the results are shown in Figures 3.6 and 
3.7 as the solid curves. The experimental runs were simulated on the basis 
of the following data. 

1. The vapor-sink driving temperature difference 

2. The numerically-integrated experimental total gas content 

3. Physical property data and conductances. 

With these input data and the governing differential equations, the gas 
concentration profiles were calculated and the gas front positions determined. 
The analytical and experimental data are in good general agreement, but the 
model rather consistently predicts a somewhat shorter gas zone. 

The general underestimation of gas zone length by the analytical model 
is not attributable to the governing differential equations but rather to a 
discrepancy between the optical absorption-derived gas density in the stagnant 
zone and the value calculated from the ideal gas law. This difference is 
summarized in Table 4.1, a comparison of heat transfer and gas concentration 
data with theory. Although the difference in concentration amounts on the 
average to only about 5%, it is the cause of the gas front shift. The heat 
transfer rate as a function of driving temperature difference, given in 
Figure 3.8 and Table 4.1, is also in agreement to within ±5%; again, the 
differences are primarily attributable to the discrepancy in gas zone length. 

In the low heat flux runs the optically-derived stagnant zone concen- 
trations are in better agreement with the ideal gas law, and the data can be 
compared with the model in greater detail, as in Figure 3.10. If axial vapor 
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TABLE 4.1. A COMPARISON OF EXPERIMENTALLY-DETERMINED VALUES 
5 FOR HEAT TRANSFER AND GAS CONCENTRATION WITH VALUES 

i PREDICTED BY THE ONE-DIMENSIONAL ANALYTICAL MODEL 


Stagnant Zone 


Sink Tempera- 

Vapor-sink 

Temperature 

Heat 

Transfer Rate (W) 

Gas Concentraj 
(std cm^/cm' 

;ion 

ture 

(°C) 

Difference (°C) 

Expt 

Model 

% Dev 

Expt 

Model 

% Dev 

14 

nom. 

35.65 

2.09 

2.05 

-1.9 

.285 

.293 

2.e 



38.8 

4.58 

4.57 

-0..2 

.330 

.340 

3.0 



41.8 

6.95 

7.00 

0.7 

.377 

.388 

2.9 



45.05 

9.37 

9.54 

1.8 

.427 

.445 

4.2 



48.2 

11.95 

12.03 

0.7 

.482 

.508 

5.4 



51.45 

14.3 

14.56 

1.8 

.565 

.579 

2.5 



58.17 

18.8 

19.24 

2.3 

.691 

.740 

7.1 

40 

nom. 

23.23 

2.22 

2.01 

-9.5 

.330 

.323 

-2.1 



26.9 

4.67 

4.57 

-2.1 

.394 

.398 

1.0 



30.8 

7.13 

7.21 

1.1 

.474 

.485 

2.3 



34.75 

9.65 

9.72 

0.7 

.567 

.582 

2.6 



38.65 

11.92 

12.33 

3.4 

.660 

.690 

4.5 



42.85 

14.3 

14.9 

4.2 

.754 

.817 

8.4 



51.1 

18.8 

20.2 

7.4 

>.993 

1.104 

• 



transport were completely neglected in the modeling, then the gas concen- 
tration would be defined by a fin- type axial temperature profile in the 
stagnant region. This profile has been superimposed on Figure 3.10 to show 
the influence of vapor diffusion into the stagnant gas. The concentration 
data of Figure 3.10 are correlated well by the one-dimensional insoluble gas 
model. Under most operating conditions, the gas concentration profiles did 
not show any gross effect of the highly soluble gas. 

There was, however, a definite "reduction" in noncondensable gas as the 
stagnant zone gas concentration increased, i.e., as the heat transfer rate 
was increased. It is very probable that this was due to gas being absorbed 
into liquid in the stagnant zone. At a sink temperature of 14°C the Ostwald 
coefficient for chlorine in carbon tetrachloride is 31.0, while it is 18.5 
at 40°C. Assuming no loss of gas from the gas zone, the total moles of gas, 
N^, must remain a constant. For a sharply defined gas zone in length, and 
for a trapped volume of liquid V. 


= (£ A 

g V 


+ aV ) C 


gvm 


(4-1) 


But in an optical absorption experiment such as done here, only the gas-phase 
noncondensable is seen, and this observable amount is 


= '<0 - Cgvm <‘•-2) 

Therefore, if stays constant, it can be predicted that the observed moles 
of gas will appear to linearly decrease with increasing stagnant gas concen- 
tration, and that the rate of decrease will be proportional to the Ostwald 
coefficient. Furthermore, in the limit of going to zero, all test runs 
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at different sink temperatures should intersect at = 0 and this will 
represent the total gas content, N^. 

Measured total contents shown in Figure 3.9 are in agreement with all 
predictions above. The linear rate at which gas is "lost" at 14°C is about 
1.9 times the rate at 40°C, and the Ostwald coefficient is about 1.7 times 
higher at 18°C than 40°C. In addition, intersection of the two linear gas 
profiles in Figure 3.9 occurs very close to zero gas concentration. There- 

3 

fore, the total molar gas content in the heat pipe is about 4.75 std cm . 

It is informative to now use (4-2) to calculate the amount of liquid 
contained as droplets and surface film within the stagnant zone. This calcu- 

3 

lation indicates a liquid load of only 0.039 cm is sufficient to produce 
the virtual gas loss shown in Figure 3.9. 

is) 

This effect was the basis of another study' ' where it was shown that 
under certain conditions a liquid absorption matrix can be used as a gas 
reservoir for gas-controlled heat pipes. This technique was shown experi- 
mentally to reduce reservoir volume by a factor of more than five compared 
to standard gas reservoir designs. 

The high solubility of chlorine in CCl^ also produced another unique 
effect not normally encountered in gas-loaded heat pipes. In the tests at 
low sink temperature, there was an unusual increase in chlorine concentration 
in the gas front region. This is shown very well in the 18.8 watt run of 
Figure 3.6. Because the heat pipe is very nearly a constant-pressure system, 
this zone must actually be cooler than the heat sink. In the fully stagnant 
zone, the wall and liquid films are close to sink temperature; a higher con- 
centration of chlorine in an intermediate zone implies a local vapor pressure 
(and film temperature) that must be below sink temperature. Another related 
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and unusual manifestation of high power level operation was observed in 

similar experiments carried out with an R-11 heat pipe and chlorine gas. In 

this case, droplets continuously formed within the vapor phase and fell 

through the gas front zone until they impinged on upward-moving vapor. This 

created a stagnation zone of high droplet density. A time-lapse photograph 

of this phenomenon is shown in Figure 4.1. The vertical streaks are falling 

droplets, while the lower diffuse disk is the droplet stagnation zone. This 

photograph was taken through the optical slit of the test apparatus in Figure 

3.2a with strong back lighting. This phenomenon is occurring within a tube 
2 

of only about 1 cm cross-sectional area. 

At the 18.8 watt heat flux for each sink temperature, the gas front was 
also very unstable with large "temporal" fluctuations in the axial concen- 
tration profile. The position of the gas front shifted axially 0.15 to 
0.20 cm in an erratic fashion with a time constant on the order of 15 to 45 
seconds, and the vapor core temperature also fluctuated about 0.1°C. 

The 14°C profile is shown in greater detail in A of Figure 4.2, along 
with a theoretical profile based on the soluble-gas modeling developed in 
Section 2.2. The axial positions of the experimental and theoretical profiles 
have been adjusted so that the profiles overlap. The experimental data range 
is indicated by the bounds given as dotted lines in A of Figure 4.2, while 
the theoretical profile is displayed as a solid line. The effects of the 
high chlorine gas solubility on the previously discussed insoluble gas con- 
centration profiles were very modest, and the differences in profile would be 
indistinguishable in A of Figure 4.2. This was primarily because of the thin 
fluid films (on the order of 0.001 cm) associated with the vertical reflux 


mode. 


FIGURE 4.1. A TIME-LAPSE PHOTOGRAPH OF DROPLET FORMATION SEEN 
AT 18 WATTS IN A GLASS HEAT PIPE CHARGED WITH Qj AND 
R-11. THE DROPLETS ARE BEING PHOTOGRAPHED WITH 
STRONG BACKLIGHTING THROUGH THE OPTICAL SLOT 
SHOWN IN FIGURE 3.2a. THE DROPLETS FALL UNTIL RESIS- 
TANCE IS ENCOUNTERED FROM UPWARD-MOVING VAPOR, 
CREATING THE DIFFUSE DISC-SHAPED AREA BELOW THE 
FALLING DROPLET STREAK PATTERNS. ONE-SECOND 
ELAPSED TIME. 





RADIAL FLUX AT INI €«F ACE GAS CONCENTRATION 
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FIGURE 4.2. CONCENTRATION PROFILES AND CALCULATED RADIAL MOLAR FLUX DENSITIES 
FOR THE 19 WATT. I 4 OC SINK EXPERIMENTAL RUN. NOTE BAND-STRUCTURE DUE 
TO FLUCTUATIONS IN GAS CONCENTRATION. 
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The unsteady gas front behavior is believed to be related to radial free 
convection. If a Grashof number is calculated based on the tube radius and 
the vapor-sink temperature difference, a value of about 170,000 is found. 

For free convection between parallel vertical plates the critical Grashof 
number is only 8000. If significant radial concentration gradients were 
present in the transitional zone from free condensation to stagnation, then 
the radial Grashof number could be high enough to initiate a convective 
motion, and this unstable convective motion would explain the concentration 
fluctuations in A of Figure 4.2. This could also explain the cold zone indi- 
cated by absorption measurements at the 14°C sink temperature. A possible 
scenario leading to this behavior will now be discussed. 

In B of Figure 4.2, the calculated radial fluxes of vapor (<(>yp) and non- 
condensable gas (<l>gR) at the liquid-vapor interface are shown. The vapor 
condensation rate shows a maximum at a point 0.2 cm into the gas zone, and 
decreases monotonically beyond this point. Chlorine gas is simultaneously 
being evolved from the film between 0 s z^ 5 0.35 cm, and dissolved into the 
film from 0.35 cm to the end of the transition zone. Both the vapor and 
chlorine radial flux densities are small but finite at 1.05 cm, the position 
of maximum chlorine concentration in A of Figure 4.2. If the vapor phase gas 
concentration rises and falls in response to the disordered mass transfer 
occurring in the first 1/2 cm of the zone, chlorine will be absorbed and 
desorbed from the liquid film. From Figure 3.4, the heat of solution of CI 2 
in CCl^ is about 14,000 Joules/g-mole, as calculated from the slope of log a 
versus 1/T, while the heat of vaporization of CCl^ is about 31,000 Joules/ 
g-mole in this temperature range. Because the film has 20 to 30 times as 
much CI 2 per unit volume as the vapor phase, it is possible that gas desorp- 
tion from the thin film under semistagnant conditions could lead to a net 
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cooling of this zone, and a relative maximum in chlorine vapor phase concen- 
tration. However, because the modeling and experimental phases of this study 
were basically one-dimensional approaches, it is not possible to substantiate 
the mechanisms described above, and these processes should be regarded as 
speculatory at this time. 

4.2 Model Validity 

The one-dimensional models developed in this study assume that axial 
transfer of energy is predominated by diffusive mass transfer and that the 
sensible heat capacity of the diffusing vapor is negligible compared to the 
latent heat capacity. These assumptions are validated in Table 4.2, which 
compares axial transport of heat via gas-phase conduction, ordinary diffusion 
and thermal diffusion. The comparison was based on the theoretical profile 

TABLE 4.2. THE RELATIVE IMPORTANCE OF VARIOUS ENERGY 
TRANSFER PROCESSES FOR A REPRESENTATIVE 
CONCENTRATION PROFILE 

Axial Heat Flux Density 


Mode 

{Vi/crt 

Conduction 

0.0073 

Thermal diffusion^ 

0.0004 

Latent heat, ordinary diffusion 

0.092 

(21 

Sensible heat, ordinary diffusion' ' 

0.015 


^Assuming a thermal diffusion ratio of 0.05. 

(21 

' 'Based on the total vapor-sink temperature difference. 
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in A of Figure 4.2. The comparison was done at the axial position with maxi- 
mum vapor temperature rate-of-change versus z^. This position was z^ * 0.38 
cm. Ordinary diffusion is shown in Table 4.2 to be the primary mode of energy 
transfer, with heat transfer by conduction accounting for less than 8% of 
the total thermal flux, and thermal diffusion accounting for less than 0.5% 
of the energy transfer. The initial assumptions pertaining to energy trans- 
port modes are therefore generally satisfied by model solutions. However, 
the energy associated with vapor phase sensible heat was neglected, yet it 
could account for up to 15% of the axial energy flux. Inclusion of sensible 
heat would be an appropriate basis for a more accurate model of gas front 
concentration profiles. 


47 


5.0 SUMMARY 

One-dimensional models for condensation in the presence of insoluble and 
soluble gases have been developed for a tubular geometry. The models have 
been compared with experimental data derived from Cl 2 /CCl^ vertical reflux 
heat pipes and are in good agreement with measured concentration profiles at 
low heat flux densities. At high heat flux densities, the gas fronts become 
unstable and exhibit a turbulent and fluctuating concentration profile. This 
turbulent behavior was associated with unusual effects, including an apparent 
depression of heat pipe wall temperature below sink temperature and the for- 
mation of liquid droplets in the vapor phase. 

The experimental concentration profiles were determined using a new 
technique having several advantages over conventional thermocouple or thermis- 
tor temperature measurement methods. By selecting a noncondensable gas which 
strongly absorbs light over discrete portions of the spectrum, the concentra- 
tion of noncondensable gas as a function of axial position was determined 
directly from optical absorption measurements in a glass heat pipe. Axial 
resolution is improved over temperature measuring methods and the vapor-phase 
flow field is undisturbed. 

Apart from the effects at high power level, gas solubility did not have 
a significant influence on gas concentration profiles. This was predicted 
by the soluble gas condensation model, and was physically attributable to the 
relatively thin fluid films (on the order of 0.001 cm) associated with the 
vertical reflux mode. The short time constants for liquid-phase diffusive 
equilibrium precluded any significant effect on gas front shape. This may 
not be true for systems with thicker fluid films or lower liquid phase dif- 
fusivities, and each case must be considered individually. 
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Liquid films and droplets within the stagnant portion of the gas-rich 
zone were found to create a virtual sink for noncondensable gas. That is, 
as vapor phase gas concentration was increased, the absorption of increasing 
amounts of gas within this liquid was detected as a 10 to 20% loss in total 
moles of noncondensable gas in the vapor phase. This phenomenon could be 
used in gas-controlled heat pipes to increase thermal control capability or 
reduce gas reservoir size by dissolving the noncondensable gas into a volu- 
metrically efficient liquid reservoir in lieu of the common gas-bulb 
reservoir. 
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APPENDIX A 

COMPREHENSIVE TABULATION OF OPTICAL 
ABSORPTION AND THERMAL DATA 


Run A1 
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Vapor core temperature: 
Sink temperature: 

Heat input: 

Calculated gas content: 
Zero-absorption current; 


49.65°C 

19.2°C 


2.087 W 
4.342 std cm^ 
4.03(x 10"^°) A 


Concentration 
3x 


position (cm) 

Current (x 10"^^)A 

(std cm^/ 

1.294 

3.960 

.0070 

1.344 

3.940 

.0091 

1.394 

3.830 

.0205 

1.444 

3.730 

.0313 

1.494 

3.580 

.0482 

1.544 

3.420 

.0674 

1.594 

3.260 

.0878 

1.644 

3.100 

.1095 

1.694 

2.920 

.1358 

1.744 

2.810 

.1529 

1.794 

2.690 

.1727 

1.844 

2.550 

.1972 

1.894 

2.410 

.2235 

1.944 

2.310 

.2436 

X.994 

2.240 

.2583 

2.044 

2.185 

.2703 

2.094 

2.150 

.2781 

2.144 

2.130 

.2827 

2.194 

2.120 

.2850 

2.244 

2.110 

.2873 

2.294 

2.110 

.2873 

2.344 

2.090 

.2920 

2.394 

2.075 

.2955 

2.494 

2.080 

.2943 

2.594 

2.080 

.2943 

2.694 

2.090 

.2920 

3.194 

2.080 

.2943 

4.194 

2.080 

.2943 

5.194 

2.090 

.2920 

6.194 

2.090 

.2920 

7.194 

2.090 

.2920 

8.194 

2.120 

.2850 

9.194 

2.120 

.2850 

10.194 

2.120 

.2850 

11.194 

2.120 

.2850 

12.194 

2.130 

.2827 

13.194 

2.120 

.2850 
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Run A2 

52.8°C 
14 \ 

4.584 W 
4.299 std cm^ 

4.03(x 10"10) A 

Concentration 

-10 3 3 

Axial position (cm) Current (x 10 )A (std cm /cm ) 


Vapor core temperature: 
Sink temperature: 

Heat input: 

Calculated gas content: 
Zero-absorption current: 


1.694 

4.080 

-.0049 

2.194 

4.080 

-.0049 

2.694 

4.070 

-.0039 

2.944 

4.070 

-.0039 

3.094 

4.030 

.0000 

3.194 

4.020 

.0010 

3.294 

3.970 

.0060 

3.344 

3.890 

.0142 

3.394 

3.780 

.0259 

3.444 

3.630 

.0425 

3.494 

3.390 

.0711 

3.544 

3.170 

.0998 

3.594 

2.910 

.1373 

3.644 

2.740 

.1643 

3.694 

2.610 

.1865 

3.744 

2.510 

.2045 

3.794 

2.380 

.2294 

3.844 

2.300 

.2457 

3.894 

2.180 

.2714 

3.944 

2.080 

.2943 

3.994 

2.020 

.3088 

4.044 

1.970 

.3212 

4.094 

1.930 

.3315 

4.144 

1.925 

.3328 

4.194 

1.890 

.3420 

4.244 

1.890 

.3420 

4.294 

1.890 

.3420 

4.344 

1.890 

.3420 

4.394 

1.890 

.3420 

4.494 

1.890 

.3420 

4.594 

1.930 

.3315 

4.694 

1.930 

.3315 

5.194 

1.940 

.3289 

6.194 

1.930 

.3315 

7.194 

1.930 

.3315 

8.194 

1.940 

.3289 

9.194 

1.940 

.3289 

10.194 

1.940 

.3289 

11.194 

1.940 

.3289 

12.194 

1.940 

.3289 

13.194 

1.935 

.3302 


Run A3 
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Vapor core temperature: 55.8°C 

Sink temperature: 14.05°C 

Heat input: 6.949 W 

3 

Calculated gas content: 4.22 std cm 

Zero-absorption current: 4.03 (x 10”^^) A 





Concentration 


[std cm^/cm^] 


1.694 

4.080 

-.0049 

2.194 

4.080 

-.0049 

2.694 

4.060 

-.0030 

3.194 

4.040 

-.0010 

3.694 

3.940 

.0091 

4.194 

4.020 

.0010 

4.694 

4.030 

.0000 

4.794 

4.030 

.0000 

4.894 

3.990 

.0040 

4.994 

3.910 

.0121 

5.044 

3.820 

.0216 

5.094 

3.630 

.0425 

5.144 

3.430 

.0662 

5.194 

3.120 

.1067 

5.244 

2.880 

.1419 

S.20A 

2.680 

.1744 

5.344 

2.540 

.1990 

5.394 

2.380 

.2294 

5.444 

2.270 

.2519 

5.494 

2.170 

.2736 

5.544 

2.065 

.2979 

5.594 

1.940 

.3289 

5.644 

1.870 

.3474 

5.694 

1.800 

.3668 

5.744 

1.780 

.3725 

5.794 

1.740 

.3842 

5.844 

1.710 

.3932 

5.894 

1.730 

.3872 

5.944 

1.730 

.3872 

5.994 

1.730 

.3872 

6.094 

1.730 

.3872 

6.194 

1.730 

.3872 

6.694 

1.760 

.3783 

7.194 

1.780 

.3725 

8.194 

1.780 

.3725 

9.194 

1.790 

.3696 

10.194 

1.780 

.3725 

11.194 

1.780 

.3725 

12.194 

1.790 

.3696 

13.194 

1.764 

.3771 




2.194 

4.080 

-.0059 

3.194 

4.020 

.0000 

4.194 

3.970 

.0050 

5.194 

3.980 

.0040 

5.694 

3.990 

.0030 

5.794 

3.940 

.0081 

5.894 

3.990 

.0030 

5.994 

3.990 

.0030 

6.094 

3.980 

.0040 

6.194 

3.970 

.0050 

6.294 

3.940 

.0081 

6.394 

3.900 

.0122 

6.444 

3.830 

.0195 

6.494 

3.770 

.0259 

6.544 

3.640 

.0403 

6.594 

3.370 

.0726 

6.644 

3.130 

.1042 

6.694 

2.820 

.1502 

6.744 

2.610 

.1853 

6.794 

2.410 

.2224 

6.844 

2.260 

.2529 

6.894 

2.120 

.2838 

6.944 

2.030 

.3051 

6.994 

1.890 

.3408 

7.044 

1.835 

.3557 

7.094 

1.730 

.3859 

7.144 

1.650 

.4105 

7.194 

1.590 

.4300 

7.244 

1.570 

.4367 

7.294 

1.560 

.4401 

7.344 

1.560 

.4401 

7.394 

1.560 

.4401 

7.694 

1.580 

.4334 

8.194 

1.600 

.4267 

9.194 

1.610 

.4234 

10.194 

1.610 

.4234 

11.194 

1.610 

.4234 

12.194 

1.610 

.4234 

13.194 

1.600 

.4267 
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Run A5 


Vapor core temperature: 
Sink temperature: 

Heat Input: 

Calculated gas content: 
Zero-absorption current: 

62.3°C 
14.1°C 
11.952 W 
4.046 std cm^ 
4.04(x 10"^°) A 

Concentration 

Axial position (cm) 

Current (x 10"^^)A 

(std cm^/cm^) 

2.194 

4.120 

-.0078 

3.194 

4.040 

.0000 

4.194 

3.980 

.0060 

5.194 

3.960 

.0080 

6.194 

3.930 

.0111 

7.194 

3.980 

.0060 

7.494 

3.990 

.0050 

7.594 

3.980 

.0060 

7.694 

3.940 

.0101 

7.794 

3.780 

.0269 

7.844 

3.640 

.0424 

7.894 

3.370 

.0747 

7.944 

3.090 

.1120 

7.994 

2.760 

.1621 

8.044 

2.475 

.2122 

8.094 

2.275 

.2521 

8.144 

2.120 

.2862 

8.294 

1.980 

.3199 

8.244 

1.873 

.3478 

8.294 

1.780 

.3738 

8.344 

1.680 

.4037 

8.394 

1.590 

.4326 

8.444 

1.500 

.4637 

8.494 

1.425 

.4915 

8.544 

1.410 

.4972 

8.594 

1.376 

.5106 

8.644 

1.367 

.5142 

8.694 

1.370 

.5130 

8.794 

1.370 

.5130 

8.894 

1.400 

.5011 

8.994 

1.405 

.4992 

9.094 

1.406 

.4988 

9.194 

1.415 

.4953 

10.194 

1.480 

.4709 

11.194 

1.450 

.4820 

12.194 

1.450 

.4820 
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Run A6 

Vapor core temperature: 65.65°C 

Sink temperature: 14.2°C 


Heat Input: 14.263 W 

3 

Calculated gas content: 3.92 std cm 


Zero-absorption current: 

3.98(x 10" A 


1 


Concentration 

Axial position (cm) 

Current (x 10"^^)A 

(std cm^/cm^) 

2.194 

3.830 

.0154 

3.194 

3.820 

.0165 

4.194 

3.840 

.0144 

5.194 

3.920 

.0061 

6.194 

3.940 

.0040 

7.194 

3.980 

.0000 

8.194 

3.970 

.0010 

8.394 

3.940 

.0040 

8.494 

3.880 

.0102 

8.594 

3.920 

.0061 

8.694 

3.930 

.0051 

8.794 

3.930 

.0051 

8.894 

3.850 

.0133 

8.944 

3.730 

.0262 

8.994 

3.580 

.0431 

9.044 

3.380 

.0671 

9.094 

2.980 

.1213 

9.144 

2.680 

.1687 

9.194 

2.360 

.2275 

9.244 

2.180 

.2654 

9.294 

1.970 

.3150 

9.344 

1.844 

.3481 

9.394 

1.740 

.3778 

9.444 

1.630 

.4117 

9.494 

1.550 

.4382 

9.544 

1.460 

.4702 

9.594 

1.400 

.4930 

9.644 

1.320 

.5253 

9.694 

1.270 

.5468 

9.744 

1.220 

.5694 

9.794 

1.210 

.5741 

9.844 

1.185 

.5859 

9.894 

1.210 

.5741 

9.944 

1.210 

.5741 

9.994 

1.215 

.5717 

10.094 

1.215 

.5717 

10.194 

1.230 

.5648 

10.694 

1.280 

. 5424 

11.194 

1.305 

.5317 

11.694 

1.120 

.5295 

12.194 

1.300 

.5338 

12.694 

1.300 

.5338 

13.194 

1.276 

.5442 

13.694 

1.230 • 

.5648 


I 
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Run A7 

Vapor core temperature: 
Sink temperature: 

Heat input: 

Calculated gas content: 
Zero-absorption current: 


72.27°C 
14. 1°C 
18.802 W 
3.534 std cm^ 
4.03(x 10"^°) A 


Axial position (cm) 

2.194 

3.194 
^•.194 

5.194 

6.194 

7.194 

8.194 

9.194 

10.194 

10.294 

10.394 

10.494 

10.594 

10.694 

10.794 

10.894 

10.994 

11.094 

11.194 

11.294 

11.394 

11.494 

11.594 

11.694 

11.794 

11.894 

11.994 

12.094 

12.194 

12.694 

13. 194 

13.694 


Current (x 10 ^^)A 

4.300 
4.100 
4.000 
3.900 
3.850 
3.820 
3.600 
4.020 
3.980 
3.980 
4.030 

4.000 
3.880 
3.780 
3.750 
3.350 
2.950 
2,500 
1.890 
1.480 

1.300 
1.160 
1.090 

.965 

.926 

.907 

.916 

.945 

.960 

1.000 

1.000 

1.000 


Concentration 
(std cm^/cm^) 

-.0256 

-.0069 
.0030 
.0132 
.0184 
.0216 
.0459 
.0010 
.0050 
.0050 
.0000 
.0030 
.0152 
.0259 
.0291 
.0761 
.1313 
.2064 
.3420 
.4696 
.5407 
.6052 
.6413 
.7121 
.7371 
. 7499 
.7438 
. 7247 
.7152 
.6910 
.6910 
.6910 


Run B1 
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Vapor core temperature: 63.35°C 

Sink temperature: 40.12°C 

Heat input: 2.221 W 

3 

Calculated gas content: 4.498 std cm 

Zero-absorption current: 4.03(xl0"^®) A 

Axial position (cm) Current (x 


Concentration 
(std cm^/cm^) 


1.294 

3.920 

.0111 

1.694 

4.000 

.0030 

1.794 

4.030 

.0000 

1.894 

4.030 

.0000 

1.994 

4.010 

.0020 

2.094 

4.020 

.0010 

2.194 

3.980 

.0050 

2.294 

3.890 

.0142 

2.344 

3.810 

.0227 

2.394 

3.660 

.0392 

2.444 

3.540 

.0530 

2.494 

3.380 

.0724 

2.544 

3.230 

.0918 

2.594 

3.030 

.1195 

2.644 

2.925 

.1350 

2.694 

2.820 

.1512 

2.744 

2.705 

.1699 

2.794 

2.620 

.1844 

2.844 

2.520 

.2022 

2.894 

2.470 

.2115 

2.944 

2.420 

.2210 

2.994 

2.370 

.2307 

3.044 

2.320 

.2408 

3.094 

2.270 

.2511 

3.144 

2.220 

.2617 

3.194 

2.180 

.2704 

3.244 

2.165 

.2737 

3.294 

2.140 

.2793 

3.344 

2.120 

.2838 

3.394 

2.080 

.2930 

3.494 

2.035 

.3037 

3.594 

1.970 

.3196 

3.694 

1.980 

.3171 

4.194 

1.970 

.3196 

5.194 

1.970 

.3196 

6.194 

1.970 

.3196 

7.194 

1.970 

.3196 

8.194 

1.940 

.3272 

9.194 

1.960 

.3221 

10.194 

1.970 

.3196 

11.194 

1.970 

.3196 

12.194 

1.990 

.3146 

13.194 

1.940 

.3272 

13.694 

1.930 

.3297 


1 


i Run B2 

Vapor core temperature: 

* Sink temperature: 

Heat input: 

Calculated gas content: 

, Zero-absorption current: 


Axial position (cm] 



Concentration 


'std cm^/cm^] 


1.694 

4.030 

.0000 

2.194 

4.060 

-.0030 

2.694 

4.060 

-.0030 

3.194 

4.060 

-.0030 

3.694 

3.970 

.0060 

4.194 

4.020 

.0010 

4.294 

4.020 

.0010 

4.394 

4.010 

.0020 

4.494 

3.990 

-.0040 

4.594 

3.980 

.0050 

4.694 

3.910 

.0121 

4.744 

3.820 

.0216 

4.794 

3.680 

.0368 

4.844 

3.530 

.0540 

4.894 

3.360 

.0747 

4.944 

3.180 

.0981 

4.944 

3.020 

.1204 

5.044 

2.740 

.1634 

5.094 

2.610 

.1853 

5.144 

2.500 

.2050 

5.194 

2.355 

.2326 

5.244 

2.265 

.2509 

5.294 

2.180 

.2690 

5.344 

2.120 

.2824 

5.394 

2.070 

.2939 

5.444 

2.025 

.3045 

5.494 

1.980 

.3155 

5.544 

1.940 

.3254 

5.594 

1.930 

.3280 

5.644 

1.890 

.3383 

5.694 

1.880 

.3409 

5.744 

1.840 

.3516 

5.794 

1.820 

.3570 

5.894 

1.800 

.3626 

5.994 

1.790 

.3654 

6.094 

1.780 

.3682 

6.194 

1.753 

.3759 

7.194 

1.730 

.3825 

8.194 

1.720 

.3854 

9.194 

1.730 

.3825 

10.194 

1.710 

.3884 

11.194 

1.710 

.3884 

12.194 

1.730 

.3825 

13.194 

1.700 

.3914 

13.694 

1.690 

.3944 


i 


Run B3 
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i 

I 

t 




Vapor core temperature: 
Sink temperature: 

Heat input: 

Calculated gas content: 
Zero-absorption current: 


70.95°C 
40.15°C 
7.13 W 

4.337 std cm^ 
4.03(xl0"^°) A 


Axial position (cm) 


Current (x 10‘^°)A 


Concentration 



2.194 

4.030 

.0000 

3.194 

4.030 

.0000 

4.194 

3.990 

.0040 

5.194 

4.020 

.0020 

6.194 

4.010 

.0020 

6.294 

4.010 

.0020 

6.394 

4.010 

.0020 

6.494 

4.020 

.0010 

6.594 

3.980 

.0050 

6.694 

3.910 

.0121 

6.744 

3.810 

.0226 

6.794 

3.670 

.0379 

6.844 

3.480 

.0599 

6.894 

3.230 

.0914 

6.944 

3.020 

.1204 

6.994 

2.750 

.1618 

7.044 

2.620 

.1836 

7.094 

2.440 

.2162 

7.144 

2.285 

.2468 

7.194 

2.120 

.2824 

7.244 

1.990 

.3130 

7.294 

1.920 

.3305 

7.344 

1.855 

.3476 

7.394 

1.800 

.3626 

7.444 

1.660 

.4035 

7.494 

1.700 

.3914 

7.544 

1.690 

.3944 

7.594 

1.650 

.4066 

7.644 

1.640 

.4097 

7.694 

1.605 

.4208 

7.744 

1.600 

.4224 

7.794 

1.570 

.4322 

7.844 

1.560 

.4355 

7.894 

1.560 

.4355 

7.994 

1.540 

.4422 

8.094 

1.540 

.4422 

8.194 

1.505 

.4541 

8.694 

1.490 

.4594 

10.194 

1.490 

.4594 

11.194 

1.470 

.4665 

12.194 

1.480 

.4629 

13.194 

1.470 

.4665 

13.694 

1.450 

.4737 


Run B4 
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Vapor core temperature: 

74.85° 


Sink temperature; 

40.1° 


Heat input: 

9.653 W 

q 


Calculated gas content: 

4.309 std cm 


Zero-absorption current: 

4.03(x 10‘^°) A 

Concentration 

Axial position (cm) 

Current (x 10 ^^)A 

(std cm^/cm^) 


2.194 

4.030 

.0000 

3.194 

4.030 

.0000 

4. 194 

3.980 

.0050 

5.194 

4.010 

.0020 

6.194 

3.990 ■ 

.0040 

7.194 

4.020 

.0010 

7.694 

4.020 • 

.0010 

7.794 

4.000 

.0030 

7.894 

3.990 

.0040 

7.994 

3.990 

.0040 

8.094 

3.980 

.0050 

8.194 

3.870 

.0163 

8.244 

3.820 

.0216 

8.294 

3„630 

.0424 

8.344 

3.400 

.0697 

8.394 

3.140 

.1036 

8.444 

2.870 

. 1428 

8.494 

2.750 

.1618 

8.544 

2.490 

.2068 

8.594 

2.295 

.2447 

8.644 

2.165 

.2723 

8.694 

2.020 

.3057 

8.744 

1.845 

.3503 

8.794 

1.740 

.3796 

8.844 

1.160 

.4192 

8.894 

1.560 

.4355 

8.944 

1.150 

.4524 

8.994 

1.460 

.4701 

9.044 

1.425 

.4829 

9.094 

1.405 

.4904 

9.144 

1.395 

.4942 

9.194 

1.370 

.5038 

9.244 

1.360 

.5078 

9.294 

1.350 

.5117 

9.344 

1.320 

.5238 

9.394 

1.320 

.5238 

9.494 

1.310 

.5279 

9.594 

1/310 

.5279 

9.694 

1.290 

.5362 

10.194 

1.270 

.5447 

10.694 

1.270 

.5447 

11.194 

1.260 

.5490 

12.194 

1.230 

.5622 

13.194 

1.220 

.5667 

13.694 

1.220 

.5667 



Run B5 
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Vapor core temperature: 
Sink temperature: 

Heat input: 

Calculated gas content: 
Zero-absorption current: 


78.87°C 
40.22°C 
11.916 W 
4.195 std cm^ 
4.03(x 10“^°) A 


Axial position (cm) 


Current (x 10 ~^^)A 


Concentration 
(std qn^/cm^) 


2.194 

4.020 

.0010 

3.194 

4.030 

.0000 

4.194 • 

3.980 

.0050 

5.194 

4.020 

.0010 

6.194 

4.010 

.0020 

7.194 

4.010 

.0020 

8.194 

3.990 

.0040 

9.194 

4.020 

.0010 

9.294 

3.970 

.0060 

9.394 

3.830 

.0205 

9.444 

3.670 

.0379 

9.494 

3.380 

.0722 

9.544 

3.120 

.1063 

9.594 

2.810 

.1522 

9.644 

2.540 

.1977 

9.694 

2.330 

.2376 

9.744 

2.120 

.2824 

9.794 

1.935 

.3267 

9.844 

1.790 

.3654 

9.894 

1.673 

.3995 

9.944 

1.560 

.4354 

9.994 

1.440 

.4773 

10.044 

1.350 

.5117 

10.094 

1.276 

.5421 

10.144 

1.230 

.5622 

10.194 

1.180 

.5850 

10.244 

1.175 

.5874 

10.294 

1.166 

.5916 

10.344 

1.140 

.6042 

10.394 

1.130 

.6091 

10.444 

1.125 

.6116 

10.494 

1.120 

.6141 

10.594 

1.110 

.6191 

10.694 

1.080 

.6345 

10.794 

1.080 

.6345 

10.894 

1.080 

.6345 

11.194 

1.070 

.6398 

12.194 

1.080 

.6345 

13.194 

1.135 

.6066 

13.694 

1.030 

.6607 


Run B6 
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Vapor core temperature: 
Sink temperature: 

Heat Input: 

Calculated gas content: 
Zero- absorption current: 


83.1°C 
40.25°C 
14.336 W 
4.157 std cm^ 
4.04(x 10"^°) A 


Axial position (cm) Current (x 10~^°)A 


Concentration 
(std cm^/cm^) 


2.194 

3.194 

4.194 

5.194 

6.194 

7.194 

8.194 

9.194 
9.694 
9.994 

10.094 

10.194 

10.294 
10.344 

10.394 
10.444 

10.494 
10.544 

10.594 
10.644 

10.694 
10.744 

10.794 
10.844 
10.894 
10.944 
10.994 
11.044 

11.094 
11.140 

11.194 
11.244 

11.294 

11.394 

11.494 

11.594 

11.694 

11.794 

12.194 

12.694 

13.194 

13.694 


4.040 

4.030 

4.030 

4.030 

4.020 

3.980 

4.000 

4.040 

4.030 
3.990 
3.980 
3.980 
3.890 
3.830 
3.690 
3.470 
3.050 
2.780 
2.420 
2.180 
1.970 
1.790 
1.605 
1.490 
1.370 
1.270 
1.165 
1.080 

1.030 
0.980 
0.965 
0.936 
0.938 
0.908 
0.895 
0.890 
0.890 
0.880 
0.880 
0.880 
0.880 
0.880 


.0000 

.0010 

.0010 

.0010 

.0020 

.0060 

.0040 

.0000 

.0010 

.0050 

.0060 

.0060 

.0152 

.0215 

.0367 

.0622 

.1172 

.1580 

.2211 

.2702 

.3191 

.3666 

.4220 

.4607 

.5051 

.5460 

.5935 

.6359 

.6620 

.6673 

.6990 

.7168 

.7155 

.7348 

.7434 

.7468 

.7468 

.7536 

.7536 

.7536 

.7536 

.7536 




2.194 

4.020 

.0010 

3.194 

4.020 

.0010 

4.194 

2.980 

.1263 

5.194 

4.020 

.0010 

6.194 

3.980 

.0050 

7.194 

3.980 

.0050 

8.194 

3.940 

.0091 

9.194 

4.070 

-.0039 

10.194 

4.050 

-.0020 

10.694 

4.030 

.0000 

11.194 

4.030 

.0000 

11.294 

4.000 

.0030 

11.394 

3.900 

.0132 

11.494 

3.900 

.0132 

11.594 

3.880 

.0152 

11.694 

3.870 

.0163 

11.794 

3.640 

.0413 

11.894 

3.300 

.0823 

11.994 

2.680 

.1733 

12.094 

2.180 

.2690 

12.194 

1.460 

.4701 

12.294 

1.180 

.5850 

12.394 

0.985 

.6859 

12.494 

0.860 

.7661 

12.594 

2.741 

.8604 

12.694 

0.650 

.9487 

12.794 

0.650 

.9487 

12.894 

0.640 

.9594 

13.194 

0.630 

.9704 

13.694 

2.610 

.9932 


APPENDIX B 


SELECTED TABULATIONS OF CALCULATED TEMPERATURE AND 
CONTENTRATION PROFILES BASED ON THE 1-DIMENSIONAL 
INSOLUBLE GAS MODEL 


NOTE: The sign appearing in the printout is a misprint for the numeral 


3". 


Run A1 
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Vapor Temperature 
Sink Temperature 
Gas Content 

Calculated Heat Transfer Rate 


49.65°C 

14.0°F 

4.342 std cm^ 
20.47 W 


Axial Position T'J’ T^ 

(cm) ■_ ! 


Gas Concentration 
(st.d ciii3/cm3) 



1 . 

1 . 0-> Of'*-, 

j . orc'.f'4 
1 . 1 0 :=■ 

1 . 1 1 
: . I '••■f r/H 


sp ciQT’ii 
r* ? '•* '"i 5 

7''" 

C | Z .' ^ ^ w . 

??=■•" ?■= 

9 1 ?7 : 


1.0 


0.0 


1.5'?9ir 


1 . ='29bf> 

. 91 194.3 

.9 95 

. 2224236— 2 

1.241 14 

. 906295 

. ?'=•= 30 3 

. ?4J 9 796-2 

1 . * 2 

. 2 9 '^4 35 

. 9 -’9 359 

. 5163796-2 

1 .224] 

..991 Oil 

. 992739 

. 7640676-2 

1 , j:55‘rif: 

. 920762 

. 97 4'92'6 

. 1 1 05*46-1 

1. 22.7 07 

. •2694 ? 

. 964419 

. 1562246-1 

1 .4lf:Sf. 

. 2599 1 7 

, 450633 

. 91 52 026-1. 

1.45 •:> 0 3 

. 9 3 b 969 

. ?3«323 

. .= 29 7756-1 

J . ii;i 51 

.917555 

. 91206:3 

. 2773356-1 

1 . 1 ? 

. 79571 

. 996969 

. 4302046-3 

1 . 544 43 

. 771559 

. 35791 1 

• T*y^*!T* 1 

1 .5759*' 

-t flr ^ 1 

. »• rr 

. 9 2.5 6 1 9 

. 72 09 266- 1 

1 . 4 0 i' 4 4 

r> 1 V 

. 79 06 0 3 

. •9525966-1 

1 . *2 3.99 3 

. 697953 

. 7=^ 36 

^ 'ic. 73^'idc - 1 

1 . r.704 1 

^ S' 

• t.* - » -• 1 1 

. 7 1 53 39 

.1121 93 

1.70: 33 

.626.272 

. b>'6791 

. 1.9524 1 

1 . 73337 

. 594924 

. 6 39255 

. 1 37'46»:. 

1 . 7*£.49.:. 

. 56 36.94 

. 600533 

. 149-4 21 

1 . 7 96 ?4 

. 529S57 

. 56596 

. 1611 21 

) . 3 27'r!2 

. 5096.97 

. 529919 

.171532 

1 . =:5-‘^9 

. 473402 

. 44527 6 

. 121157 

1 . 99 079 

.445125 

. 4^3439 

. 190022 

i. • r'w c w 1 

.41797 

.4 34 .: 37 

. 19917 

1 . - f. 3 75 

. 39200? 

.404263 

. =05651 

1 . 9 25 .2 3 

. 367275 

. 379025 

. 21 2517 

?. 01-7? 

. '•4.3796 

. 2 52 '9 01 

. 

. * • .* 2 c' 

. :- .215 26 

. 32 9-:4 3 

. 224^ 05 

0794 9 

. 300504 

. 207053 

. 2294 : 

11117 

. 2 ‘2 066 

.296211 

. 2349 06 
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Run A1 (continued) 



.361967 

. 26669 

. 335? 

r . 1 74 1 3 

.£44331 

. £43404 

. £4 :-i4 35 

i . c’ 056 1 

.££7357 

. £«1£95 

. 347343 

i?. s' 371 

. £1£?47 

.£15367 

. 350751 

•='.?655i 

. 1 37301 

. £00 £*36 

. £539£4 

.£.30006 

. 1£417£ 

. 136333 

. 35t966 

I. 33154 

.171411 

. 1 73193 

. £59715 

£ . 3 6 5 03 

. 15947 

. 160379 

. 362257 

.£. 3 34f.l 

. 14.3305 

. 149573 

. £646 01 


. 137371 

. 1 33939 

. C‘*iO 1 t o 

-.45747 

. K=-31£5 

. 123 001 

. £68767 

£ . 43336 

. 1 1 90=7 

. 1 19745 

.37 06 15 

£. 53044 

. 1 1 05 33 

.111113 

.37£3£4 

3.55133 

. 10£6£ 

. 1030 = 3 

' .373303 

£.5334 

. 95 £3*333:- 1 

. 955'='73£-l 

.£75363 

•3. 61433 

.333605E-1 

.956337E-1 

. £76714 

3. 64637 

..319536E-1 

.831457E-1 

. £77963 

3.67735 

. 75938 lE-1 

.761133E-1 

. £79117 

:•. 7n ;» .:!6 

. 704957E-1 

. 705046E-1 

. £*? 0 1 £i6 

.3.74 0S3 

.65£699£-l 

.653901E-1 

.£81174 

3. 77c 3 

. 6 04654E— 1 

.604443E-1 

. £92088 

3.? 0373 

. 559936E- 1 

.5594£7E-1 

.. £82934 

•3 . :? 3536 

.51S473E-1 

. 5i?6££E-3 

. £8 £7 16 

£. 96675 

.47S3C5E-1 

.47331£E-1 

. 3944 S*5| 

c . 9 93 £3 

.444035E-1 

.443794E-1 

.£951 08 

3.93971 

. 04 1 093 

. 040933 

.£85727 

?. 461 19 

. 379966E- 1 

.?79?41E-1 

. £86299 

£. 9936*3 

. 351331E-1 

. 034966 

. £86? £8 

?. 0£416 

.3£4773E-1 

. 3£«fi3!E-l 

.287317 

3. 05564 

. 30C14TE-1 

. 39335*E-l 

.237769 

■ . HrTII; 

. £773 £3 '• -1 

. 37550=6-1 

. c-28 166 

1 1 3r 1 

.£561 746- 

. 354'? *496- 1 

. ~ 

1 =*'irr< 

. ? -6534'=- ! 

. £34 7-" _'6-l 

. c 933 r'3 

■-t. 1 -•! 

. £*1 S445*---l 

. Sit 6336- ■ 

. S^’^35- 

. c 1 Or 

.£0l6‘=':-t£-1 

. 133-7-6-1 

. 36 35'-: 

1*4 4 '“■4 

. 1 9* 1 j 1 

. 01'=4-7 

. 399-1--* 

£ . 76 

.171 73..-6-1 

. 1 7ti'.'"66- 1 

. 33i'l O'-' 

:-.075 

. 01 6 .-: 44 

. 1 5r. 7:1 = ► - 1 

, S* :<ri* 41; 

=1. :-5£'-=' 

. 1 4t 147E-1 

. 1^4-.4£*=-l 

. £'“05--.’ 

3704 7 

. 1 34r?-4---i 

. :3'-.£4?6-T 

. I30‘’i7 

- .401 9*" 

. 1 J 

. 1 333 0- 6- 1 

. 3=* ••9“’ 5 

■*. 4:-;7a? 

. 1 J 45366- : 

. 1131t-»6-1 

. 33 1 1 3 :• 

3. 4643? 

. 1 056 - '-‘6 — 1 

. 01 Oa = 

. £313=- 

•* . 4 4 

. £735? 7f- £ 

. 06«E-= 

. £3H £f 

3. '^•3 7 S3 

. ?37.£59F-r 

. S3-4 3S46 — 3 

.3315. 7.= 

? . 55'^ 

. 9 : 6 =• - £ 

.3161S16-? 

. £9 16 36 

v.59'tS5 


. 7507S16-? 

. £-=131 : 

..6e'-'33 

. 701944f-£ 

. 691 473- -£ 

. £9191 ^ 


I 

t 
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Run A1 (continued) 



■ 6' 4 6 7 ) *6 F — 3 

. 63*- 935E-3 

.393017 

'• . t-Sf'c '■* 

.5957-^3E-3 

. 5-?,Ai)9?P-? 

.393107 

? . r i t r i 

. 5*4 3'374 £-3 

. 54 012 9E--' 

. 3991 

74 - Cl 

. 50565.'-f:-r' 

. 497443£-3 

. .69 3367 

•>. ?y?74 

. 4r:.5957fc-9 

. 459164F.-3 

. .39 3 3 37 

:■< . n i£‘ 

. 4 393 1 H f-3 

. 4.33 01 7E-: 

. •3994 1'i 

>•. -.4?7J 

. 3?-=-4:£{9^-:‘ 

1 ^ c cr — I 

• •» r r 1 T ’T ^ 

. 3^:46) 

. =! 74 1 

. 36’*4 5 ?fc - 3 

. 35 1 5“ — 3 

. e. ^ "'T' 1 6 

V. . •f'Otif.7 

. :j:-:591 3 £-3 

. > ?i'07 ?f-3 


?. ?371? 

. 3n’967-:F-3 

. 304337E-3 

. .:"? 36 1 

** • 7 

.^?5557t-3-_ 

.3905-3^-3 

. 39t*:’54 

4. nodi ? 


. 3597066-: 

. 3 '9 36 93 

4. t i}:l4 

.34-3053E-? 

. 3 «9 69 96 -3 

. •393739 

4 , 06 ?• 06 

, 3343793:-:' 

. 330 344 ?,—? 

. •3‘*-'3763 

4 . (•••4457 

. 3n7349t:-3 

. 3035156-3 

. 393793 

4 . 1 4,£. rif, 

. 19154SF-3 

. 1 99 O'? 3£— 5 

. 39391 =• 

4 . 1 5753 

. 1771969-3 

. 173979E-3 

. 393344 

4. 1 ,?'4ni 

■ 164 0fj6E— 5 

.161 0556-3 

. 3"9 3'367 

4 . 3305 

. 15r?7=^E-3 

. 1 4934-9E-3 

. 3 9 3? 3 9 

4 . 35 1 9‘z' 

. 14C993-3 

. 1 3 04736- 5 

. 393903 

4 . 3334^ 

. 1309756-c 

. 133656-3 

. 393935 

4. ■:1495 

. 131 96? -3 

. 1 1 971 36-3 

. 393941 

4. ?4643 

. 1 13591 £-3 

. 1 115966-3 

. 3 '93956 

4. 3779) 

. 1 06033-3 

. 1 043436-3 

. 393969 

4 . 4 0339 

, ?imoii4c -5 

. 9760 36- 3 

. 3 93'9'? 

4 . 44033 

. 993053E-? 

.9163766-3 

. 3'93991 

4.47336 

. 9774 9 :-?-:-! 

. 36 374 9.6-3 

. 393001 

4 . 5|* ;'94 

. 9?975f.-3 

. 915063--;' 

. 393 1! 09 

4 . 5 35 33 


. 7739946 -? 

. 3931(17 

4 . 5t69] 

. 7479"lE-3 

. 7«593?£-3 

. 3 '9"'!0c'3 

4 , 5 '4'4.5 ';( 

. 7 1 535 3 

. 70??0*:-6-3 

. c"?30?9 

4. 3:977 

.68737E-? 

. ‘-('-•suit — :- 

.39? ''34 

4 . "tic 1 35 

. 1 ii4£- •; 

. : 535 I 36-« 

. 39.30 :9 

4. -9374 

. :45939£-3 

. 6 3 5 3'9?i6 - 3 

. 39? 041 

4. 73433 

. 63079;E- 3 

.6313E-5 

. ?9:m;'44 

4. 7557 

. 6305 093-"“ 

.6 1 1 13 36-3 

. 39«'.i45 

4. 79719 

. 6 1 4 . 3694 - > 

. 6051316-3 

. 39 3 ''4 7 

4.91367 

.619.997E-3 

.609131E-3 

. 39 9047 
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Run A3 







■-! n ^ 

. 1 ' 

■f 1 

. ■4 4 r 

* » y 1* *r* c 

. 

5? 

• I' '-I 

^5?Sl 

. f* •' 1 C’ 4 ^ 


.? 13' ^ 

'V ? * ’’ ■? 

. 44r;7 

(' ?!H.r f 


I'r 1 74 

. '^0'3C>? 'd 
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1 

. •? ?5I:8 

1 
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1 - ~*r.r 


1 H 1 -•• 

.84:£-0; 4 

r f u*: 4 

CO O ^ 1 cr 
• »T*C C • A -• 

w r-c * 1 L 

. :::fir.404 


. r75:>04 

J4 y. i'* 0 

. 748 £6r 

;• ^ r 

. 7l-=Ti = 7 

4 I ! -• Of- 

.808887’ 

4 


4 7 1 0,? 

. 8 c. 38 1 r 

'-•I'l?':- 

. fi 8! 1 89 
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f.rf*4K 

. •58877:8 

? r**s 

.4*5708 

*7 - C ^ *4 

. 4^.f,f;79 

o f' ^ 1 

. 4;5b03 

^ '1'-* 

. 4 098f'4 


. «9?004 

7^' 4 ;l»r« 

. ?f-81 01 

r?5ii4 

. 3«4f.4 

:‘\73I 

. 3183CJ 

-:4-f!r:l 

.891?48 

•! •. f ' "• 
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?l 1 7? 

. 8f8180 
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RfyROT^UCTBILlTY OF Til’: 
DAr,R IS P(X)R 


i"C 

i°C 

std cm^ 
W 


Gas Concentration 
(std an3/ciw3) 



■U ’p» S 

. ■:07084 


. 0 0^98 - 

987f-“:-7 

. 76074? 

981 084 

. 115741 

97178-1 

. 171 47=’ 

959 1 0 9 

. 84 €798 

9 484b- ■ :' 

. 3 444 '-<8 

?81 398 

. 4c-59 0 =* 

89573 

. b-1 04 4 7 

Z- ^*TO .'•c' 

. 775877 

831 '•■8-8 

. >•=^51 1 *- 

794 -.b 1 

. 1 1 46 hk. 

75-4 -'17 

. 1 ?41 f 

714)r.:8 

. 15 34 

875? cc 

. 17 80 8 

b- 5 ^ 1. T/ 
Tt‘ •’ C .' . * i. 

. 1 =-?75-. 

5 -p p 0 p -p 

. 806861 

5551b- 

. ; -8 1 5 1 

519 004, 

. 2?54‘4d 

4-4 74.C: 

. 849864 

458479 

. 8599 =•? 

4821 1« 

. 8704 86 

8 9 •' b- H 7 

. 88 0186 

868877 

. 899909 

3415; 5 

. 286^19 

319:388 

. 304233 

898 4 8 4 

. 5 1 0-4 1 9 

875908 

.'*1 7 0 35 

*-56 708 

. 5? =‘6 39 
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Run A3 (continued) 
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Run A3 (continued) 
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. 190 354E-8 
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. 17E.i:)1E-.? 
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. 337565 
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. 1 ri4 4'?c E— c‘ 

. 39768 1 

147ms 

. 

. 975^.65E-? 
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1 77 s,’ 
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. ? r'4c7£E- !•! 

. S?7’t.75 
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. 6-1 78£- - 3 
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. 5-'=‘144E-3 
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. 5 9 77 £5 
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. 3977£‘- 
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Vapor Temperature 
Sink Temperature 
Gas Content 

Calculated Heat Transfer Rate 


63.35°C 
40.12°C 
4.49 std cm^ 
2.014 W 


Axial Position 

t;;- 


T 

f 

(cm) 

w 




1 . 

• ?3”65 


1. 

0 

] . T - : -4 1 

. ? 3f' 




1 . 3 ' t'-i4 

. '4 -‘ 4? 4 




J . 7? -4 7 

.?34:?i 




i . 7r 

. =• ??*:• 




: . 7=*0^’.3 

. ?c'377ir. 




I . 0'=- 

. 3 i ~ 




1 . 

. ??0.?'’- 




1 . ? r -■ 1 1 

.=^16147? 




1 . ? 1 4 

. 9161 MI 




1 . J 7 

. 91 3-='9r 


. 9-5 

1 . •4737 

. ?m5-77 


m ^ m 

■'c:965 

r‘. Ofi4E ? 

. (14 5 I I 1 


. r--03lT. 

3. 1134 76 

. '=. 9?4 =5 


• 

•• 655 :' 

3. 045 ?■? 

. r 9 ? 0 0 1 


. 9 

■: 1 7 4 

c'. 0 ^ 5.33 

. ?65I-77 


. ?757 i-' 

3 . 1 £v -•4 

.97941 1 


. 9 (V 769.3 

f. i^'-:?7 

. 5 5 f 9 7 


. 9= 

' r 1 

?. 1 -74 

. 4? 


. 44 

y ff- 1' r* i’ 

3.317-='? 

. ?4 0 : 41 s 


.i " 

■3?1 1 

: . 34 

. 9649?4 


. 91 

579 

?. r,7?? ? 

.•:n767= 


. 5'“ 

* *£ 7 ^ 

•• . : 09z<c 

. 7 r* r s ? 


. r 7 

5375 

i m ^4 1 • 05 

. 766-69 



1 434 

?. 770=^7 

. 74443 4 


. 33 

5415 

3. 4111 1 

• 7 c’ 1 c' 4 


. 7 9 

7534 

3 . 4 ? 1 6 '? 

. t 9 069 9 


. 1 ' 

0079 

?. 463 !f 

. t-715 719 


. 9 9 

“44 :■ 

3 . 4 ? .£6r* 

.647343 


. 71 . 


3 . T'?^3i3 

. 63 L'6 64 


.67 

4 09 

3. •fo?75 

■ *S »“i 


. t.4 

t Oof. 

c.f-3 46? 
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. 6 I 
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?. r 14? 1 
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. 57 

•?94 

3. -.45^4 

. 5 1 ?9 0 ? 


. -?4 

9379 

3. 67566 
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. 5! 

r 4- 7 

3. 7 m- JSC* 

. 4635.57 


• 4 r» 

9‘r”4 

3 .7 71 ; 6 :• 

. ;!9 m :-! 1 


• 

1 6 7? 

. 75745 

. 41 4 1 : f. 


• 4 

5399 

c . 7--7=*=- 

.3911 = 


. 4 0 

•■^77 


Gas Concentration 
(std cin3/cm3) 

0.0 


. :•=•! 1 n jH-* 

. ar.df. 

• T' r.' I ^ ^ 

. 1 ij H4 -.4 

. I 3Sc-r.^? 

. in 


- 1 

-1 
_ 1 


. ru 1 r 


. 5 'll 1 -1 

. 
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. 1 1 j 3 
. 1 34=.3?, 
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. 14 7=*^! 

. 1 7 7 0 3 - 
. I 3T-f;73 
. 1 =*f.riT-» 

. 3'!??7? 

. *• 1 t i 
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Run B1 (continued) 
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.36?£43 

1 o??r 

. 31 0374 
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Run B1 (continued) 
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4. r:7?95 

.99979J9-9 

. 90555:-F-9 

. -9934 

4 . 4044? 

.779054E-9 

. 99.7'?? 9E-? 
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. 175493E-9 
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. 999957 
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. 15! ?4 3E-9 
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9 Or* 74 

. ^.47994P- 9 

. 144431F-9 
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5> . 94 0 9 9 
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. 1 994.99E-? 

. 3 £ 3“ 03 

?. croPc 

. 1 z' 4-4 l.i4 — — c! 
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-'Ol 95 
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. 1 99 7E4E-9 
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- . 9t.?4J 

. 1 99997E- 9 

.1971 01 E- 9 
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. 1 94545E-9 
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. 1 9595 • E-9 
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Run B2 


Vapor Temperature 
Sink Temperature 
Gas Content 

Calculated Heat Transfer 


Axial Position 
(cm) 

T* 
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^ . ? r> r 

. ‘43 

-1 . ? '•' 5 4 

. 45^3y;- 

4. 

. 4u'4”rj 

4. 

. 4 1 i f 

4 . 41 --C .4 

. ’433?: : 
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4. f f t f = 

. :* i 1 1 f :■ 

^ 1 T' * 


4 . 4 j 11 
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.•4 1^-:--’ 

4.4.?.: If 

.'411 jit: 

• r* * ^ ^ ? 

. •?f.ri':.4 

4 . ?.:-:4i-=- 

. '? 0 r oi' £' 

*♦ . r 7 ? 

. r 

4 . ?'4’=;;T 

. z'9'S :^- ': 

4 , c 5 5 

. '1 ?‘r*“^i:i5 

4 . •ff.'f ?{4 

• r f 4 (1 5 i 

■i. 

. 350 4 5 1 

4 . 1 

. 4.ir:?r. ;-. 

4. 4 47 4.5 

. r i f 4.^ n 

4 , 'r* ? ■• ■ 4 

. iCfnOf 


. rS4:-iS 


.7f7rif 

!• “ r<iia 

.744 = .45 

1 f"i?? 

.7l4Pi 

. 1 1 1 

. f4aoi3 

f j f n 

.►-.713 

1?415 

. f. 

• 5 . i ?c’f .ri 

.31144^ 

f.. ?5331 

.533571 

5 . .E:3.??4 

.555T.17 

^ ) 4 c.' ? 

.547745 

j44t< 

. 5 01*1445 

7* — . ^ -• 

.4 7.?7f £• 


. 4 47??'4 

^ y-z 34 

.4 ?574t 

4r.f.'f‘l 

■ . * 335 71 

, *4 ^ 4 <4 

. :753f 7 

1.' e? .• 

. T c . -j <' 

. 3531 f 

•T C“ C C' C 
• • - * * ■ •. 

.331 -?e€. 
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Run B2 (continued) 
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. £' ? S 4 9 9 
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77 ?£ 

.s 11 0 :- 1 
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('ir 7 ? 
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. t T * 1 r 

. 17£?77 

. 1 7571 9 
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. i-:.o^94 

. 1 8 594 

. 54 0.? 9 

7 . ? r ^ :• 

. 1 ?. 0 ? f 1 
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. 55 15 5:? 

f> . 

. i 4 0 .*;t.? 
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. r'S'5 ?l 

. i c0c.4f- 

. l?55-5 

. 35 5 1.55 

r . 0 1 ^ >+ 

. lac'Cit. 

. 155791 

. 38 1 04.-: 

o. 'i4*f r*-. 

. 1 1 S7:-:r 
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-.. or 74^ 

. 1 05'=*?S 

. 1 1■•7.::55 

. 5 8 8 1 7 4 

t.. 1 Oc'O-i 

. ??7?1 ?E-1 

. 999.857E-1 

. 38.54 77 

*:.• . 1 ,* ri^'5 

. 91 =>4 0iE-l 

. 95951 45 - 1 

. 570819 

. 1 r. ? 0 -- 

. 0?7'7'9c: 

. S’r.4'^57E- i 

. 3758 1 3 

r.. 1 - -7 1 

. 79E4I.I??- J 

. 804:: 1 

. 37448'^ 

t.. rVO 14 

.74C'^92E-1 

. 74775SF-1 

. 5781 95 

o. .?r.i>:.7 

. 4.89 1 94.E- 1 

. 854907E-1 

. « 7 79 OH 

■r . c.’ ■? 1 c 

. C6407 

. 845855E-1 

. 379595 
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. r-?54 9:-:E-1 

. 599855E-J 

. ::?0iS.94 



.558795E-1 

. 531978 

•f*. 

.5159‘=-.0£-l 

. 518S:48E-1 

. 395 178 

''.41 r! ? 1 

. 4 7724.5E- 1 

. 47985'iE- 1 

. 5 ’5 4:? 9 5 

T- . 4 4 -•‘:'.4 

. 44:: 081 E- 1 

. 444 =*55E-1 

. ”9555 

4. 47457 

.411J89E-1 

.41 -'889E~ 1 

. ?!98595 

7-04 ? 

. 03315 

. or-! 5.58 5 

. 55719 

r'.. t. 4 '5 

. 353854E-1 

. •54r.59E-l 

. 3 5 5055 

4. 

. 3331 19E-1 

. 035585 

. 359795 

*T • ^ *r^ *z* 

.“0417EE-1 

. t;04485E-l 

. 1P59515 

7 . 7. £-7 01 

.E"319 9-:.E-1 

. 5‘519.5'8E-1 

. 39'<! 75 

- .6?7T'4 

. 0981 IS 

. 581 095E-1 

. 5907 =»-. 

4. ft? 07 

. 941 ■=•:■ VP- 1 

. 54188=E-1 

. 391389 

- . 7 1 f r 

. 394097E-1 

• 5 53 885F— 1 

, ■-'■Q 1 .j Cl 1 

4 . 74 4 t 

. c 074*;p8-i 

. 5 (i8 ^ ? 1 E ~ 1 

. 599 5 '-^4 

- . 7 7 ^ r. 

. t r*194.-E-l 

.1^14 "5t- ! 

. 3 5.5555 

. 7 1 < ' 1 4 

. 1 7789 VP- } 

. 1 77001E-: 

. 59 357- 

•■ . ■*'4h7 ,:> 

. IC.4 537E-1 

, 18 584V8-J 

. 59j--.-.? 

:■ ..471 -• 4 

. 151E58E-1 

. 15] 557E-1 

. .vvjrtjr? 

•: . ?01 77 

. O 1 4 Cif.a 

. 1 i'9 77 9'=.-l 

. 594?r •* 

. 4 : 5 - 

. J 1 

. 1591 4.^p-l 

. 59J8- 1 

-. . ■ 

. 1900-‘4E-1 

. 1 I 5 5 r 7 E - J 

. 3 5 4 =* 7 


• 1 1 09.:‘cE- 1 

. 1 1 017^8-1 

. 555 .: 3^ 

0::‘54 r 

. 1 r954?E-l 

. 1 01 759E-1 


7. 0^.4 4c 

. 009475 

. ^?^‘l 7F-8 

. 35571 5 



77 



?. 0S4<5f: 

.■374 4 4E-3 

. 3^ 71 Of 6-=- 

. 33“'-=' 54 

7. 1 154S 

. 3 03041 E- 5 

. •?00?4 7E-? 

. 336 1 1 4 

7. 14*E.0J 

. 74o£53F.-3 

. 733643E— 3 

. 3463'43 

7. 1765? 

.6330:3E-£ 

.601 653E-3 

• 3-46466 

7. 50706 

. 6361 76E-5 

.65^05E-3 

. 3366 1 4 

X • ■- •• 1 T* 7 

.53735SE-3 

. 530436E-3 

. ?-467 41 

7 . 56 ? 1 5 

.545c36E-'' 

. 5«r56l3E-3 

. «36833 

7. 54565 

.500633E-3 

. 434573E-5 

. 347013 

7.3541? 

.46335-3E-3 

. 45617E-3 

. 337134 

7.35471 

. 43636 lE-3 

.451 066r —3 

. 537536 

7.3405^ 

. 3-44377E“? 

. 338746E-3 

. 53733x 

7. -IS 077 

. 364376E-3 

. 353006E-3 

.53740- 

7.451S4 

. ?;.‘67746-5 

. 531 665E-? 

. 5374?^ 

7.4313c- 

.311 33'^E-S 

. ?0653-?t-5 

. 33756 

7.51335 

.3S8033E-3 

. 583473E-3 

. 337637 

7.5433? 

. 36670-4E-3 

.36333E-3 

. 347663 

7.57?41 

.3471 llE-3 

. 343-461 E-5 

. 337746 

7.60334 

. ;3-4175E-=- 

. 335844E-3 

. 337737 

7.63447 

. 31i?SbE-3 

. 3 03 Or- 4E -5 

. 537844 

7.665 

. 137333E-3 

. 1 4’i 3 1 46—3 

. 537887 

7.635‘=-.? 

. 1 34c36€-c 

. r503='7E-£ 

. 347456 

7 . 736 05 

.171 33-=*6- r. 

. 1T.3735E-3 

. ?=*7363 

7.75653 

. 1607136-5 

. 157717c-£ 

. 347494 

7.7-371 1 

. 150643E-3 

. 1477-43E-3 

. 54805 5 

7. ?J76-a 

. 1416144.-3 

.138303E-3 

. 336043 

7.64317 

. 1 33556E-3 

. l«03?lE-3 

. 34807= 

7.8737 

. 136413E-3 

.13344^6-3 ' 

. 338035 

7. 4 0-433 

. 150157F-3 

. 117784E-5 

. =*431 1 

7.43376 

. 1 14736E-5 

. 113443E-3 

. 3481 36 

7.-57033 

. 1 10033E-3 

. 107833E-C 

. 338153 

3. 00031 

. 1 06££:E-5 

. 1 0404E-3 

. 33815 

4. 031 34 

. 103033E-3 

.10101 3E-C' 

. «381 53 

06137 

. 1 00673E-3 

. 3364 1 7E-3 

. 3'46l 66 

6. C-334 

.-4S3663E-3 

. 4A‘:iS'44C’- 3 

. 3381 71 

5. 13533 

. 973436E-3 

, 4S-iS46F- 4 

. 338174 

?. 15346 

.376034E-3 

. 3563 05€- 3 

. 3381 75 
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APPENDIX C 

LISTING OF COMPUTER CODE FOR 1-DIMENSIONAL 
MODEL OF CONDENSATION HEAT TRANSFER IN THE 
PRESENCE OF AN INSOLUBLE GAS 


NOTE: The sign appearing in the printout is a misprint for the numeral "3". 
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• ; 0-&-PEP-7b 

• M » ?b4 • * C t D ‘ £b4 • t P ‘ 5(s4> 

.-n^l < T) »i=i4*EyP ♦T*«h J 

PNP0<T» ■'il *EyP'Hj/T) ♦7**i^3 


REPRODUClClI^rJi'Y tM;’ 'u 
ORIGINAL PAGE W P0'J»5 


41' 

-iTrt 


ZO 

1 ^. 


! TaTi4i. CCNt’ENSEP LENI?TH 


. J 
2 b 


ftPQP. 

P n ■4TM0tPMEPES) 

IN FPEE 2. 


■•41' *31 .•52.>41 t!.12 !(5-PILf1»G-riN> 

. 418«»« . 033?;. . 0117* .4«ti;'? 

PEFiIi Vt.V2.V3 !i4-ViHPOP.CDf<jT. IN C*ri.HT. QF 
r"4Tt4 1 . 047.8. 81E-3. 31 000. 

35 PF^iti h0.82.i 43 !l.N<Pi» HO ♦ He-'T ♦ H«*lN<T> 

36 DHTH 30.33.'’5.-4544.47.-2. 301 34 

40 PEHD N.22.M !N»i;H'-2aNF. I N'‘.c . . 22-1' IN. LENGTH. N» IN'! P. 

41 rii=iTfl 114.. 8 .?. to. 

45 P2H E5 SNHX. EPPOP PEPMITTEp [N T— .iHu. PEP ITEPATIiDN 

4'. ["4ri4 . 0"1 
51 INPuT“T-VHPaP»"! T 1 
INP'JT"T-! IN'- »•• S T3 
lNPi.lT"5Tr'. 'ZC. 3 F '3i43»"S80 
60 H»Z? n; Nt «rt+ 1 : mi i : M3*M*''<+ 1 

35 •41*273. lS*EyP >Hi'» : P1«=NP0‘ T 1 ♦273. i . i H4 *h t 'Pi 

70 MO*T >♦ 173 . 1 51 i.l’SmT j - T 3: C5«'J4 4. : i‘ ; *H«»2*'31 ♦'.1 l*i.'2 ' • 1 ♦'’2*V? > 

71 2 1 * i *Z'P .'.11 '52 '*'?5*1 .'^..'*'31 + 1. ''j2'1'43*'^3 *jc 

72 tl *-2. - « 1 . ♦'51 ' '52 ' *‘4** ?! tc»-‘?l '52 

~3 H0*21*H 

4fi "Imt 3n < ' ll ■ 

?t H'4T i.ls . i.'js I *..1 

4-*. V ' 1 •». 4'fi*.!.: ■ 1 •»ii. : fu 1 >*V ' 1 • !>4r'.''.'rT fop CHlC. 3p C'2>.1"2> 

100 •^OP 1*2 TO N? 

10*. V ' I ' *F 'iP 1 — 1 •♦H' jt.'.I ♦>! * ■'.' 3 • V • I ' ! 'lE T I 


53 

53 


110 P • ? • * I . H 

115 FOP 1*3 TO M? •p-mhTPI ' ftp '.imH. TENP. ;'3L'.'TIC»< 

l?n p. I , ' f PI _p . i-n . j .(EX'T 1 

150 T3«o. J .5P TO 280 
I'sn r3»o. 

200 FOP 1*2 in Ml ICHLC. ■: HmIi l. ■'IhTPIX-: FQF .yHPnt TEMP. :OLN. 

205 T *'.'0^'jP* J . I I s C E ■'42 Tx ■h4*F '"P * L ♦ T* • 'Z 7*H 3— 'ZE- 

210 ; *C X* I V • 1 ♦ 1 I -./ ■ I - 1 > * ♦ • *Z 7 • 1 . - ' • ♦ ZE 'Z 7- . 5 • ' r 

il5 FmiZF*! Z'P ■ T» • < 1 . -1.' 'Z7 

230 I*- *^2. ♦Fi - • 1 . -E •E'l-l* 

235 :• I •■• 1 .♦£» 'f: t" I •■-•►♦’." I ♦ '•■♦' t . -E» ♦D' I-l ' ' P 
240 NEXT I 

250 V • Ml . ■ . F*'.i I m 3» ♦P. •t'.'l' • • ■ 2. ♦F^F. *'Z ' N ' • 

2E0 FCP I*H TO 2 ZTEP -1 ! v’PP'^f TEMPEFhTL'PE ppOPIL? 

2-?5 V' n *0' I > -c ‘ I • *v • 1 ♦! • : M= r i 

280 ri 2>*>t'2-i43. PI 

8.35 FOF I«‘: to H IP MHIPP' FOP '.ill TEMPEP HT'.'PE CHLCZ. 

230 Im 1 > ■ 'P 2 -I" I-t > » *F. I > sneXT 1 
iJ'jl.' FOP I»Mt TO M2 

30*. p . I > m I pc'*V < I-M> -p. 1 - 1 I )*F f I i : ..PWT I 

510 iJ > m3> ■ • I2*V'M1 > -2. ♦!' • M2> ' ' • B 1 - 2. pF 'M2* » x m4»'4 2 

525 FflP I ■M2 TO 2 5TEP-1 

330 T I *'4 * 1 > s 14 ^ 1 1 *P . I » -P I ) . 1 ♦ 1 > : T2*Mp3 ' T 1 -*4 • I ‘ ■ 

035 IP T2 T3 '50 TO ;4'. 

340 T3*T2 
•'4 5 'IrPT I 

400 IF T3 E5 GO TO 1 =•" 

800 r 2m <p 1 -Ff(pO <1.(0 • * (.MIS Tml IM♦"4•\/ > 1 I S 'Z • 1 > m 1 P 1 -PtiPO ' T ■ > T 
80^ I ml IO*'.'P*V '*(1 'X'" •Ml'm'Pl -♦ NP'O < T • * T 
810 5* . 5* ' Z ' 1 * ♦'! ' Ml • > X M£m . 5« 1 (,( • 1 I *14 iM 3 > > 

E12 FOP 1*2 TO Ml ! I.NrF.JPHTe N'. G HMp TEMP. PPOFIlE? 

-.13 '32* '2^I<" 1 ' X MF T I 

■bl5 FHP !■» TO M ! juTFGPHTP 5HZ H,*P TPMP. PPOFI'.e: 

8? I. T ■'.(('♦'J'-m'V *T'X Z*I'*'P1 -FMP'O . T • » T 

-.25 *5*'^'Z I'X iZ'2*'.'2^'.' • 1 ♦M ' X NC; T I 

8Z<o 5*Vl*'4"*i’X P',2*'.''4*G2*H0*(^2x 27* 'G0-.5* • < vl ••Z 3 ' 

b 4 0 75*Z*0- • r7^22*n ♦H'.*M 1 1 r.i mu«<#.3c;mi.( . | * #75 x f* Z *■*• 1 ♦.,'•2 

E4*. POP 1*1 TO M3 

©*■•0 p.I'*-r^HO» i-i :'«EXT I 

670 PPINTxPPINT 'rp.)' cNTPhTIDm HNP TEMFfpmT'.ipE ppoF ILE Z " x PP InT 
E71 ^PtHTipPlMT 'M.'3,-.(3 ?jPPINT 

8?0 Frjp i.t TO .1 

840 PPIMT p' I • .'4. I ■ IMfXT 1 

©41 PriP ImMl TO M3 

838 P'lMT p. I > I i . V I-*1' .'Z • I-M. x'lEXT I 

700 FOP 1*1 TO Ml 

710 P' I ■ *M0* < t- 1 1 : *(EPT I 

720 PPIMTxFP I'iT 

730 POP 1*1 TO Ml 

740 ■X'PIHf < I I .i4iM^I » . V< I * . : • I > IMfiT I 
3'30 MHT PPINT P*N1». 

400 CMP 


